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Abstract. Precise proton and neutron form factor measurements at Jefferson Lab, using spin observables, 
have recently made a significant contribution to the unraveling of the internal structure of the nucleon. 
Accurate experimental measurements of the nucleon form factors are a test-bed for understanding how the 
nucleon’s static properties and dynamical behavior emerge from QCD, the theory of the strong interactions 
between quarks. There has been enormous theoretical progress, since the publication of the Jefferson Lab 
proton form factor ratio data, aiming at reevaluating the picture of the nucleon. We will review the 
experimental and theoretical developments in this field and discuss the outlook for the future. 
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1 Introduction 


E2 
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One of the fundamental goals of nuclear physics is to un¬ 
derstand the structure and behavior of strongly interact¬ 
ing matter in terms of its basic constituents, quarks and 
gluons. Quantum chromodynamics (QCD) is the theory 
of the strong interaction, responsible for binding quarks 
through the exchange of gluons to form hadrons (baryons 
and mesons). The electromagnetic form factors are among 
the most basic quantities containing information about 
the internal structure of the proton and neutron, together 
known as nucleons. The challenge of understanding the 
nucleon’s structure and dynamics has occupied a central 
place in nuclear physics. High energy electron scattering 
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provides one of the most powerful tools to investigate the 
structure of nucleons. 

Early electron scattering experiments with nuclei were 
motivated by a need to verify predictions of the then cur¬ 
rent models of the electromagnetic interaction of electrons 
with nuclei, and in particular with the proton and neutron; 
Rosenbluth predicted that high energy electrons would be 
scattered dominantly by the magnetic moment of the pro¬ 
ton [1|. Available accelerators in the early fifties had en¬ 
ergies smaller than 50 MeV, and provided information on 
the nuclear radius of elements from Be to Pb. The first 
clear evidence that the proton has a structure was ob¬ 
tained at the High Energy Physics Laboratory (HEPL) at 
Stanford in the period form 1953 to 1956, under the lead¬ 
ership of Robert Hofstadter [2]. A proton charge radius 
of 0.77 fm was extracted by Chambers and Hofstadter [3] 
from the electron-proton data obtained using the electron 
beams with energies up to 550 MeV at the HEPL, confirm¬ 
ing that the proton has a finite size. In his review paper 
Hofstadter [5] discussed in detail the extraction of proton 
charge and magnetization radii between 0.72 and 0.80 fm 
using different models. Almost sixty years after the work 
of Hofstadter [J, the question of whether the proton ra¬ 
dius is 0.8775 (51) fm, the CODATA [5] value from ep 
elastic scattering, or is 0.84087 (39) fm [6], the muonic 
hydrogen value, which is a difference of seven standard 
deviations, is being discussed intensely. 

A similar change in accepted concepts occurred when 
the data, which was obtained at the Thomas Jefferson 
National Accelerator Facility (TJNAF) or Jefferson Lab 
(JLab) for the proton’s electric to magnetic form factor ra¬ 
tio, Gep/Gmp from double polarization experiments and 
completed in 2000 at a four momentum transfer squared, 
Q 2 , of up to 5.6 GeV 2 [7ll8ll9l[T0] . differed drastically from 
the form factor results obtained with the cross sections 
data using the Rosenbluth separation method Emma 


contributions to the cross sections. A direct measurement 
of the two photon contribution to the double polarization 
observables was undertaken at JLab. The results confirm 
the expectation that two photon exchange affects the pro¬ 
ton form factor ratio at less than the percent level. Com¬ 
bining these various experiments and theoretical calcula¬ 
tions, the possible role of two photon exchange in bridging 
the gap in the extraction of the proton form factors from 
cross section and double polarization observables can be 
investigated. 


1.1 History of Elastic Electron Scattering on the 
Nucleon 

Elastic electron proton scattering has evolved since the 
history making series of experiments with electron beams 
of the HEPL at Stanford in 1950s. Under the leadership 
of R. Hofstadter, a series of crucial results were obtained 
from cross section measurements [123] . Several fundamen¬ 
tal pieces of information were established following these 
experiments, including the approximate 1/Q S decrease of 
the cross section with Q 2 , establishing the approximate 
shape of the charge distribution, and a first value for the 
proton radius. Theoretical work evolved in parallel with 
these experimental “firsts”, leading to the description of 
the elastic electron scattering in terms of the lowest order 
process, the exchange of a single virtual photon with nega¬ 
tive invariant mass squared; this lowest order contribution, 
also called the Born term, was expected to be dominant 
because of the smallness of the electromagnetic coupling- 
constant oiem- Fundamental expressions for the hadronic 
current and the definition of two invariant form factors, F\ 
and F 2 , later named the Dirac and Pauli form factors, of 
the Born term, were issues of this period. In 1957 Yennie, 
Levy and Ravenhall 1241 derived an expression for the ep 
cross section in terms of these two form factors, F\ and 


14.15,16,17,18,19,20,21,22 . The standard form factor datable following Rosenbluth’s work Q], as given in Eq.[3] 


up to 1990’s had been entirely defined by cross section 
measurements, and suggested that, for Q 2 < 6 GeV 2 , the 
ratio ^ p Gep/Gmp ~ 1, where p p is the magnetic moment 
of the proton. The double polarization experiments at 
JLab demonstrated that p- p Ge p /Gmp decreased approxi¬ 
mately linearly with Q 2 for Q 2 > 0.5 GeV 2 , dropping to 
a value of 0.35 at Q 2 = 5.6 GeV 2 which was the highest 
Q 2 investigated at that time. 

In the last several years, the field of nucleon struc¬ 
ture has been investigated extensively. Two recent exper¬ 
iments at JLab have increased the Q 2 range of the Geu 
and p p Ge p /Gmp data. These data have triggered much 
activity in the determination of the flavor separated form 
factors of the dressed up- and down quarks in the nucleon. 
Experiments at Mainz, JLab and MIT-Bates have resulted 
in a much better coverage of the low Q 2 range for Ge p , 
contributing to the intensive discussion of a possible dis¬ 
agreement between electron elastic scattering and muonic 
hydrogen data used to determine the proton radius. Sev¬ 
eral experiments have measured the e + /e _ cross section 
ratio for the proton with a level of precision that tests 
recent calculations which included two photon exchange 


The possibility of measuring either the polarization 
transferred to the recoil proton, or the asymmetry if the 
target proton or neutron is polarized, with longitudinally 
polarized electrons, was discussed in a paper by Akhiezer 
et al. [25] in 1957. It was to be more than 30 years be¬ 
fore such experiments, which require a polarized electron 
beam, could be performed with good accuracy. Further pa¬ 
pers on double polarization experiments followed, includ¬ 
ing Scofield [26], Akhiezer and Rekalo [271128] . Dombey 
[25] , and Arnold, Carlson and Gross [50] . 

The construction of the Continuous Electron Beam Ac¬ 
celerator Facility (CEBAF) at JLab, in Virginia, led to an 
intensive program of nucleon form factor measurements, 
first for the proton and then for the neutron, and a sig¬ 
nificant breakthrough in our understanding of the proton 
structure. In the Born approximation the transferred po¬ 
larization has only two non-zero components, both in the 
reaction plane defined by the beam and scattered electron, 
one along the recoil proton momentum, and the other per¬ 
pendicular to it. For the proton, polarization transfer has 
been used most often at JLab mmm-, it requires a re¬ 
scattering of the proton to measure its polarization. For 
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the neutron, target asymmetry has now been used success¬ 
fully j31ll32lf33] for the determination of the Geu/Gmu 
ratio; this requires a polarized target of either 2 H or 3 He, 
which limits the maximum electron current that can be 
tolerated without significant depolarization of the target. 

Among the earliest electron-proton scattering polariza¬ 
tion experiments is a search in 1963 at the Orsay linear ac¬ 
celerator, for a one-photon/two-photon interference effect 
with an un-polarized electron beam and an un-polarized 
target. They searched for normal and transverse polar¬ 
ization components. The normal polarization component 
was found to be 0.040 ± 0.027; the transverse polarization 
component was 0.000 ± 0.028 at Q 2 of 0.61 GeV 2 [31] . 

A similar single-spin experiment in 1970 with an un- 
polarized 15-18 GeV electron beam at the Stanford linear 
accelerator, and a polarized proton target with polariza¬ 
tion perpendicular to the reaction plane to characterize 
the interference of the two-photon exchange with the sin¬ 
gle photon exchange (Born) process, produced asymme¬ 
tries of order 1 to 2 % in the range of Q 2 0.38 to 0.98 
GeV 2 [IS]. 

On the neutron side, the pioneering experiment of Madey 
et al. performed the first recoil polarization measurement 
of GEn at a Q 2 -value of 0.255 GeV 2 in 1994 [36] at the 
MIT BATES Linear Accelerator; and the first double po¬ 
larization measurement of the proton form factor ratio, 
Ge p /Gm P i was also carried out at the same lab in 1994- 
1995 by measuring the two polarization transfer compo¬ 
nents Pf and P t , at Q 2 -values of 0.38 to 0.50 GeV 2 [37] 

IMiM!. 

Also the mid nineties saw a number of double polariza¬ 
tion experiments at Nationaal Instituut voor Kernfysica 
en Hoge Energie Fysica (NIKHEF) [30] and Mainz Mi- 
crotron (MAMI) piTlPTS] to determine the neutron electric 
form factor up to Q 2 ~1 GeV 2 . All experiments used po¬ 
larized electron beams and a polarized target, either 2 H 
or 3 He. 

In this review, we will focus on the space-like nucleon 
form factors, as in the past 15 years they have been stud¬ 
ied more extensively both experimentally and theoreti¬ 
cally, compared to the time-like nucleon form factors |43j . 
Also the strangeness form factors are not discussed in this 
review; see Ref. 33] for a review of the field of parity vi¬ 
olating electron scattering and strangeness form factors. 

This review is organized as follows. Section [2] describes 
the formalism of elastic electron scattering on the nucleon. 
In subsection nm the use of elastic differential cross sec¬ 
tion data to extract the two electromagnetic form factors 
of proton and neutron by the Rosenbluth, or longitudi¬ 
nal and transverse (LT)-separation technique is reviewed. 
Subsection [U discusses how the form factors are mea¬ 
sured in the double polarization experiments. The two 
photon exchange formalism is explained in subsection l2.3l 

Section[3]is devoted to discuss the experimental status. 
Subsection 13.11 describes the experiments which extracted 
the electric and magnetic form factors for the proton and 
the neutron from measurements of cross sections. Sub¬ 
section 13.21 reviews the status of the double polarization 
experiments and discuss the results for the proton and the 


neutron, obtained from the recoil polarization method and 
beam-target asymmetry measurements. The role of two 
photon exchange contributions in the elastic ep reaction 
in resolving the discrepancy between the proton’s elec¬ 
tric form factor extracted by recoil polarization versus the 
Rosenbluth separation technique are discussed in subsec¬ 
tion [T3] Subsection 13.41 reviews the present status of the 
proton charge radius. Subsection 13.51 discusses the present 
status of flavor separation of nucleon form factors. 

Section E] deals with the theoretical interpretations of 
the electromagnetic nucleon form factors. Subsection 14.11 
reviews the models of the nucleon form factors. These 
models include conformal fits to the form factors, vector 
meson dominance, dispersion analysis, constituent quark 
models, pion cloud models, transverse densities, and cor¬ 
respondences with higher dimensional theories. Subsec¬ 
tion EE2 describes the Dyson-Schwinger equations and di¬ 
quark models. Subsection l4.3l discusses links between deep- 
inelastic scattering and nucleon form factors which in¬ 
cludes perturbative QCD inspired models and general¬ 
ized parton distribution (GPD) models. Subsection 14.41 
describes lattice QCD calculations of nucleon form fac¬ 
tors. 

Section [5] summarizes the current issues and challenges 
in the area of electromagnetic form factors. This section 
closes with a discussion of future experiments at Jeffer¬ 
son Lab which will measure the proton and neutron form 
factors to Q 2 =10 GeV 2 or greater. 

2 Formalism of Elastic Electron Nucleon 
Scattering 

The lowest order approximation for electron nucleon scat¬ 
tering is the single virtual photon exchange process, or 
Born term. The Born approximation is expected to pro¬ 
vide a good lowest order description of elastic eN scatter¬ 
ing (with N = p,n) because of the weak electro-magnetic 
coupling of the photon with the charge and the magnetic 
moment of the nucleon. The amplitude for the process is 
the product of the four-component leptonic and hadronic 
currents, and J t ,. and can be written as: 

iM = 

= l9 ^ V [ieu{k')Yu(k)\ \-iev{p')r lJ ‘{p',p)v{p)} ,( 1 ) 

where fc, k',p,p' are the the four-momenta of the incident 
and scattered, electron and proton, respectively, con¬ 
tains all information of the nucleon structure, and g is 
the metric tensor. To insure relativistic invariance and the 
correct parity property of the amplitude A4, can only 
contain p, p’ and 7 ", besides numbers, masses and Q 2 , 
defined as, Q 2 = — (q 2 — w 2 ) = — q 2 , is the negative of 
the square of the invariant mass, q of the virtual photon 
exchanged in the one-photon approximation of eN scat¬ 
tering. 
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The most general form for the hadronic current for the 
spin 4-nucleon, satisfying relativistic invariance and cur¬ 
rent conservation, and including an internal structure is: 


The Rosenbluth separation technique takes advantage 
of the linear dependence in e, in the reduced cross section 
a re d, based on Eq. ©, as follows: 


J M 


iev(p') 


(Q 2 ) 


ia^q v 
2 M 


k 3 F 2 (Q 2 ) 


v (p), 


( 2 ) 


where M is the nucleon mass; Kj , with j = p,n is the 
anomalous magnetic moment, in units of the nuclear mag¬ 
neton, hn = eh/(2M p ). The Dirac and Pauli form fac¬ 
tors, Fi(Q 2 ) and F 2 (Q 2 ) are the only structure functions 
allowed in the Born term by relativistic invariance. As is 
now the most frequently used notation, KjF 2 with j = p,n 
will be written as F 2p and F 2n , respectively. In the static 
limit, Q 2 = 0, F\ p = 1, F 2p = k p = 1.7928 and F\ n = 0 
and F 2n = n n = —1.9130, for the proton and neutron, 
respectively. 


2.1 Cross Section Experiments 


The Lab frame differential cross section for detection of 
the electron in elastic ep or en scattering is then: 


da 

/ da N 

) ^ ( 

df2 e 


Mott Fbeam ' 


(f 2 (Q 2 ) 


+ r 


F 2 (Q 2 ) + 2[Fi(Q 2 ) + F 2 {Q 2 )] 2 tan 2 y]),(3) 


with r = Q 2 /AM 2 . The Mott cross section is: 


da 

df2 


a 2 cos 2 | 

Mott 4-E/ eam sin 2 


(4) 


The incident electron (beam) and scattered electron 
energies are labeled Eb eam and E e , respectively. The frac¬ 
tion Ee/Ebeam in Eq. |3]) is the recoil correction to the 
Mott cross section. 

Experimental cross section data are most easily ana¬ 
lyzed in terms of another set of form factors, the Sachs 
form factors Ge and Gm mm- The relations between 
Ge and Gm and F\ and F 2 for proton and neutron are: 


F^E(p,n) F[ (p.n) ^~F2(p,n) 

Gm(p,ti) F\(jj n ^ T F 2 (^ p n p (5) 

The scattering cross section Eq. © can then be written 
in a simpler form, without an interference term, leading 
to a separation method for G 2 E and G 2 M known as Rosen¬ 
bluth (or Longitudinal-Transverse) technique, as will be 
seen below. Now the cross section is: 


e(l + t) Ebeam ( da_\ . ( da \ 

ared ~ r E e \dn)J / \dn) Mott 

= G~ m + ~G e , (8) 

showing that a re d is expected to have a linear dependence 
on e, with the slope proportional to G 2 E and the intercept 
equal to G 2 M . 


2.2 Double Polarization Experiments 

In 1968 and 1974 Akhiezer and Rekalo E2HE1 discussed 
the interest of measuring an interference term of the form 
GeGm by observing the transverse component of the re¬ 
coiling proton polarization in ep —> ep at large Q 2 , to ob¬ 
tain Ge in the presence of a dominating Gm- In a review 
paper Dombey m emphasized the virtues of measure¬ 
ments with a polarized lepton beam on a polarized tar¬ 
get to obtain polarization observables. Also later in 1981 
Arnold, Carlson and Gross m discussed in detail, that 
the best way to measure the neutron and proton form 
factors would be to use the 2 H(e, e'n)p and 1 H(e, e'p) re¬ 
actions, respectively. 

Indeed, both the recoil polarization and target asym¬ 
metry measurement methods have been used successfully 
to measure the proton and neutron form factors to high 
four momentum transfer, Q 2 , at JLab. The same methods 
have been used also at MIT-Bates, MAMI, and NIKHEF, 
to make precise proton and neutron form factor mea¬ 
surements at lower Q 2 . Both methods are discussed be¬ 
low, with benefits and drawbacks of using polarized target 
and/or focal plane polarimeter. 


2.2.1 Recoil Polarization Method 

With a longitudinally polarized electron beam and an un¬ 
polarized target, the polarization of the incoming elec¬ 
tron is transferred to the nucleon (proton or neutron) via 
exchange of a single virtual photon as shown in Fig. [1] 
For elastic ep scattering, in the single photon exchange 
approximation, with a longitudinally polarized electron 
beam, the only non-zero polarization transfer components 
are the longitudinal and transverse, Pi and Pt- The nor¬ 
mal polarization transfer component, P n , is zero. For sin¬ 
gle photon exchange, the transferred polarization compo¬ 
nents can be written in terms of the Sachs form factors 
as: 


da 

df2 e 


da 

df2 


Ee 


1 


Mott 


Ebeam 1 H - T 


(g 2 e + -Gm 


( 6 ) 


where e is the polarization of the virtual photon defined 
as: 


1 

1 + 2(1 + t) tan 2 ^ 


(7) 


I 0 Pn = 0 

t d 7 7-> {Ebeam H” E e ) / — ■ r 2 ^2 

IoPe = hP e - — -v t(1 + r) tan y G M 

I 0 P t = —hP e 2\J r(l +t) tanyG_EG M 

lo = G 2 e H— G 2 m 


e = 


(9) 
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Fig. 1. Illustration of the kinematics and polarization of the 
recoil nucleon for ep —» e'p and en —> e'n. 


where h = ±1 are the beam lielicity states, P e is the mag¬ 
nitude of polarization, and 9 e is the electron scattering 
angle. 

The ratio of Ge to Gm is then directly obtained from 
the ratio of the two polarization components I\ and Pi 
as: 


Ge _ Pt (Ebeam + Pe) 9 e 
~Gm ~ ~~P(. 2 M an T' 


( 10 ) 


The double-polarization approach to obtain the ratio 
at the large momentum transfer by measuring two polar¬ 
ization components simultaneously was first proposed at 
JLab in 1989 [47]; it is based on a combination of spin 
precession in a magnetic spectrometer and using a proton 
polarimeter. The major advantage of the method, com¬ 
pared to cross section measurements, is that in the Born 
approximation, for each Q 2 , a single measurement of the 
azimuthal angular distribution of the proton scattered in 
a secondary target gives both the longitudinal, Pi, and 
transverse, Pt, polarization. Thus the ratio of electric to 
magnetic form factors of the proton is obtained directly 
from a simultaneous measurement of the two recoil po¬ 
larization components. The knowledge of the beam polar¬ 
ization and of the analyzing power of the polarimeter is 
not needed to extract the ratio, Ge/Gm > strongly decreas¬ 
ing the systematic uncertainties. The kinematic factors in 
Eq. (HOD are typically known to a precision far greater than 
the statistical precision of the recoil polarization compo¬ 
nents. 


2.2.2 Asymmetry with Polarized Targets 

It was discussed by Dombey [29] in a review paper in 
1969 that the nucleon form factors can be extracted from 
the scattering of longitudinally polarized electrons off a 
polarized nucleon target. In the one photon exchange ap¬ 
proximation, the elastic electron nucleon scattering cross 
section can be written as a sum of two parts: S, which 
corresponds to the unpolarized elastic differential cross 
section given by Eq. ©, and a polarized part, A, which 
is non-zero only if the electron beam is longitudinal po¬ 
larized mm ; 


(ii) 



Fig. 2. Illustration of the kinematics and orientation of the 
target polarization P, for the reaction en —¥ e'n and ep —» e'p. 


The polarized part of the cross section, A, with two 
terms related to the directions of the target polarization, 
P{9* (/)*), is given by [151135] : 


A = —2 


/ da_\ E e 9e_ j t 

\df2 J Mott Ebeam 2 \ 1 + T 


^ sin 9* cos c/)*GeGm 


^Jt[1 + (1 + t) tan 2 y] cos 9*G 2 M S j 


( 12 ) 


where 9* and <f>* are the polar and azimuthal laboratory 
angles of the target polarization vector with q in the 2 
direction and y normal to the electron scattering plane, 
as shown in Figure [2] 

The physical asymmetry A is then defined as 


A- a +- = - 

a+ + a_ S’ 


(13) 


where cr_|_ and cr_ are the cross sections for the two beam 
helicity states. 

For a longitudinally polarized beam and polarized tar¬ 
get, the measured asymmetry, A meas , is related to the 
physical asymmetry, A, by 


Ameas — 9,P e PtargetA, (14) 

where P e and Ptarget are electron beam and target polar¬ 
ization, respectively, and A can be obtained using Eqs. © 
and (fl2l) . 


A = - 


2\J t( 1 + r) tan 


G% 


T S-l2 


sin0* cos (/)*GeGm 


r[l + (1 + r) tan 2 -d-] cos 9*G 2 M 


(15) 


From Eq. (USD: it is apparent that to extract Ge, the 
target polarization in the laboratory frame must be per¬ 
pendicular with respect to the momentum transfer vec¬ 
tor q and within the reaction plane, with 9* = 7t/2 and 


ah = E + hP e A. 
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<fi* = 0° or 180°. For these conditions, the physical asym¬ 
metry A in Eq. CSD simplifies to: 


P, = ^r~e ^ 

&red 



A 


perp — 


—2-^/r(l + t) tan 

(£e -) 2 + £ 

' Gm ' + e 


§_e Ge 
2 Gm 


(16) 


As ( Ge/Gm) 2 is quite small, A perp is approximately pro¬ 
portional to Ge/Gm- In practice, the second term in Eq. (fT5l) 
is not strictly zero due to the finite acceptance of the detec¬ 
tors, but these effects are small and depend on kinematics 
only in first order and can be corrected for, so the ratio 
Ge/Gm is not affected directly. One can also note that 
A pe rp ~ Pt / ( h l 'fi ) ■ 

The discussion above is only applicable to a free electron- 
nucleon scattering. For a quasi-elastic electron scattering 
from a nuclear targets, like 2 H or J He, corrections are re¬ 
quired for several nuclear effects. 


+ 2G m ^(SGm + A—^F 3 ) 
1 + £ M z 


(23) 


While the Sachs form factors depend only on Q 2 , in the 
general case the amplitudes depend also on e. The reduced 
cross section and the polarization transfer components, P t 
and P( are sensitive only to the real part of the two-photon 
amplitudes. The normal polarization transfer component, 
P ni is sensitive to the imaginary parts of the two-photon 
amplitudes. In the Born approximation, only the first term 
remains from a re d , Pt and Pe in Eqs. ED, (G2D and (l23l) 
while P n is zero. 


3 Experimental Status 


2.3 Two-photon exchange 

In the one-photon exchange process, the form factors de¬ 
pend only on Q 2 but not on other kinematic variables. A 
deviation in the form factors from constant when varying 
the kinematics (i.e. Eb eam and the scattered electron an¬ 
gle, keeping Q 2 constant) would indicate the presence of 
a mechanism beyond the Born approximation. 

In the general case, elastic eN scattering can be de¬ 
scribed by three complex amplitudes mmm- gm, ge, 
and i* 3 , the first two chosen as generalizations of the Sachs 
electric and magnetic form factors, Ge and Gm, and the 
last one, F 3 , vanishing in case of Born approximation. The 
reduced cross section, a re d can be written as: 


The structure of the nucleons has been investigated ex¬ 
perimentally with rigor over last 70 years using elastic 
electron scattering. The two Sachs form factors, Ge and 
Gm, required to describe the nucleon charge- and magne¬ 
tization distribution have been traditionally obtained by 
cross section measurements. In the static limit, the pro¬ 
ton, Gmp, and neutron, Gmu, magnetic form factors are 
equal to the proton, ji p . and neutron, (i n , magnetic mo¬ 
ments while the proton, Ge p , and neutron, Geu, electric 
form factors are equal to unity and zero, respectively. The 
earliest experiments at low Q 2 found: 


E , G Mp 

Gr Ep ~ 


Gmu 


(24) 


O’red — G 2 m + ~G\ + 2 GM^- {&Gm + 

+ 2 e -a E s ( sg e + f 3 ), 


where 


(17) 


The Q 2 dependence of these form factors can be approxi¬ 
mately characterized by a dipole form factor: 



0.71 


-2 


(25) 


and 


m M {Q 2 ,e) = G m {Q 2 ) + mG M (Q 2 ,e) 
MG E (Q 2 ,e) =G e (Q 2 ) + mG E (Q 2 ,s) . 


V 

AP 


s — u 
4 M 2 


\/r(l 


t) 


1+6 
1 — £ 


(18) 

(19) 


( 20 ) 


The polarization transfer components can be written as: 


Pr, = 


2e(l + e) hP e 


Gred 


- G m Z(6Ge + j^F 3 ) 
2 £ 


+ G e Z(6G m + ——F 3 )\ (21) 


Pt = - 


2e(l - e) hP e 


GeGm 

Gred 

+ G E $t(6G M ) + G m ^(6G e + ^F 3 ) 


( 22 ) 


The data for Gm p have shown good consistency between 
different experiments up to 30 GeV 2 , however, the deter¬ 
mination of Ge p at Q 2 greater than 2 GeV 2 has suffered 
from large error bars. The neutron electric form factor, 
Geu, is small and difficult to extract from cross section 
experiments. New experimental methods using spin ob¬ 
servables were needed which pushed the development of 
polarized targets and new accelerators with high duty fac¬ 
tor and polarized electron beams. 

The recent generation of electron accelerators with high 
polarization and high current electron beams, at MIT- 
Bates, MAMI and JLab, have made it possible to investi¬ 
gate the internal structure of the nucleon with precision. 
In particular, the new series of experiments that measured 
spin observables, like beam-target asymmetry and recoil 
polarization, have allowed experiments to obtain the pro¬ 
ton and neutron electromagnetic form factors accurately 
to large Q 2 . In this section we describe the proton and neu¬ 
tron form factors obtained from cross sections and double 
polarization experiments. 
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3.1 Cross Section Experiments 

The electric and magnetic form factors of the nucleon can 
be extracted from measurements of cross sections at a 
constant Q 2 and different beam energies using Eq. ©. 
This is known as the Rosenbluth separation technique. 
The strict linearity of the reduced cross section is based 
on the dominance of one-photon exchange in the elastic 
electron-nucleon scattering reaction. In principle, only two 
e points are needed to determine the slope, G 2 e /t 1 and 
intercept, G 2 M , from Eq. ©. Usually experiments have 
measured more than two e points for a given Q 2 . From a 
practical experimental viewpoint, more e points allow bet¬ 
ter understanding and checks of systematic errors. From a 
theoretical viewpoint, the linearity of the e dependence of 
the reduced cross section can be investigated with more e 
points. One clear sign of a two-photon exchange contribu¬ 
tion to the cross section would be a non-linearity in the 
e-dependence of the reduced cross section. Unfortunately, 
the two-photon exchange contribution can have a linear 
e-dependence which cannot be experimentally separated 
out in the cross section measurement and must be calcu¬ 
lated theoretically. The form factor data presented in the 
following sections were not corrected for hard two-photon 
exchange contributions when they were extracted from the 
elastic scattering measurements. 

3.1.1 Proton Form Factors 

Extraction of the proton form factors from the cross sec¬ 
tion data is complicated by the strong dependence of the 
Mott cross section on the scattering angle. In addition, 
as can be seen in Eq. ©, the relative contribution of the 
two form factors to the reduced cross section changes with 
Q 2 . This difficulty in measuring both form factors within 
the same experiment with small errors bars was a problem 
from the first experiments using the Rosenbluth separa¬ 
tion technique. Ref. m contains a tabulation of measure¬ 
ments of Gmp and Ge p from the early 1960’s which are 
plotted as open triangles in Figs. [3] and 0] From Ref. [IT . 
we have selected to plot only data which have relative error 
bars of less than 10 %. One can see that Ge p is measured 
with this precision only to Q 2 <0.2 GeV 2 while the Gm p 
data points are plotted only above Q 2 = 0.2 GeV 2 . 

Throughout the 1960’s, experiments increased their 
precision. Measurements were done at the Stanford Mark 
III accelerator of cross sections to 2% statistical preci¬ 
sion at Q 2 between 0.18 to 0.8 GeV 2 [TO] . Their Ge p and 
Gm p values are plotted as multiplication sign symbol in 
Figs. [3] and |U In 1971, results on Ge p and Gm p were pub¬ 
lished which combined cross section measurements at the 
Cambridge Electron Accelerator with previous cross sec¬ 
tion measurements to cover a range of Q 2 between 0.13 to 
1.75 GeV 2 . The results are plotted in Figs. [3] and 0 as open 
circles [T3]. There is good agreement with the Ge p and 
Gm p of Ref. [T21 i 11 the region of overlap, while Ge p /Gd 
for Q 2 > 1 GeV 2 starts to drop-off below unity. In 1970, a 
SLAC experiment measured Ge p and Gm p at Q 2 between 
1 to 3.8 GeV 2 and found Ge p /Gd > 1, which is opposite 


to the trend of Ref. m- These data are plotted as filled 
diamonds in Figs. [3] and [4] In this same time period, Ge p 
and Gm p measurements were done at Bonn [T5 for Q~ 
between 0.34 to 1.94 GeV 2 and at DESY ljl for Q 2 be¬ 
tween 0.67 to 3.0 GeV 2 . The data from Ref. X5] and [TO] 
are plotted in Figs. [3] and 0 as filled square and crossed 
diamond, respectively. Both of these data sets agree with 
the downward trend in Ge p /Gd for Q 2 > 1 GeV 2 seen 
in Ref. [T5]. which disagrees with the rise in Ge p /Gd ob¬ 
served in Ref. [T4j . 

In the 1970’s, a series of experiments at Mainz sought 
to measure the Ge p and Gm p form factors at Q 2 below 
0.1 GeV 2 with greater precision. The first El measured 
cross sections at Q 2 between 0.014 to 0.12 GeV 2 . The 
second |T8i did measurements up to Q 2 of 0.055 GeV 2 with 
an emphasis on extracting the charge radius of the proton. 

The Ge p and Gm p form factors of Ref. \T7\ and ifS are 
plotted as crossed square and open square in Figs. [3] and 0 
These experiments demonstrate the precision that can be 
obtained in the measurement of Ge p at extremely low Q 2 . 

In Sec. 03 results from these experiments for the proton 
charge radius are shown in Fig. 1161 

In the late 1960’s and early 1970’s, experiments at 
SLAC [53 pushed the limits of the ep cross section mea¬ 
surements to Q 2 = 25 GeV 2 . In the 1990’s, an experiment 
[541 at SLAC measured the ep cross section at forward 
angles for Q 2 from 3 to 30 GeV 2 with improved statistical 
precision. Both experiments extracted Gm p under the as¬ 
sumption that pl p Ge p /Gm p = 1- The data from Ref. [53] 

(open squares) and [53 (open stars) are plotted in Fig. 0 
The data are consistent with each other and show a drop¬ 
off in Gm p /p p Gd above Q 2 = 7 GeV 2 . 

In the 1990’s, at SLAC, two experiments were done 
which extended the precision and upper range of Q 2 for 
measurement of Ge p by the Rosenbluth separation tech¬ 
nique. Ref. [[TO] measured at Q 2 of 1.0, 2.0 and 3.0 GeV 2 . 

Ref. [20] measured at Q 2 between 1.75 and 9 GeV 2 and 
this experiment pushed the measurement of Ge p to the 
maximum Q 2 that has been done at this time. In Figs. 0 
and 0 the data of Ref. [TO] and [20] are plotted as filled 
stars and open diamonds, respectively. The Gm p data of 
both experiments agree with each other and previous ex¬ 
periments, while the Ge p /Gd are very different at Q 2 = 3 GeV 2 . 
The Ge p /Gd of Ref. [30] agrees well with the early mea¬ 
surements of Ref. [13; at Q 2 = 1.75 GeV 2 and have a flat 
Q 2 dependence with a slight rise for Q 2 > 1.75 GeV 2 . A 
global reanalysis of cross section experiments was done by 
Ref. [52] in 2003. This reanalysis found that the cross sec¬ 
tion measurements from different experiments were con¬ 
sistent with each other in extraction of both Ge and Gm, 
though it excluded the small angle data of Ref. [TO] from 
the global analysis. 

The Ge p /Gm p data measured by the recoil polariza¬ 
tion method will be discussed in the upcoming Sec. 13.21 
but the effect on the extraction of Gm p will be briefly dis¬ 
cussed here. In 2002, an extraction of Gm p was done from 
the cross section data of the previous experiments using 
the constraint that 

Ge p IGm p = TO - 0.13 x (Q 2 - 0.04), 
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which originates from a fit to Ge p /Gm p from the recoil po¬ 
larization experiments [56]. The fit to the extracted Gm p is 
shown in Fig. 0] as a dashed line and the fit is a few percent 
larger than Gm p from the standard Rosenbluth method. 
Also shown in Fig. [I] as a solid line is a fit by Ref. [55] 
to the Gm p extracted by Ref. [56 with additional low Q 2 
Gm p values from Ref. i58j. Since the time of these fits, the 
effects of two-photon exchange on the extraction of Ge p 
and Gm p have been calculated by theorists and this topic 
will be discussed in Sec. 13.31 

Taking advantage of the high duty factor of modern 
accelerators, cross-sections can be measured to high pre¬ 
cision over a range of e in a relatively short time period. 
This was recently done with spectacular precision by an 
experiment at MAMI. Cross sections were measured at 
1422 kinematic settings covering a Q 2 range from 0.004 to 
1.0 GeV 2 with average point-to-point systematic error of 
0.37% [52]• Data were taken with the three separate spec¬ 
trometers of MAMI at 6 different beam energies. With 
this large data set, the authors extracted Ge p and Gm p 
by fits to their cross section data rather than the tra¬ 
ditional Rosenbluth separation technique. The group has 
published a long paper [50] on the same data set and in¬ 
cluded other world data in their fits. The sensitivity to 
different functional forms for the fits was investigated by 
using many different spline and polynomial forms. The 



Q 2 (GeV 2 ) 


Fig. 3. Gep/Gd extracted from cross section measurements 
versus Q 2 . The data from before 1980 are: open triangle (red) 
ITT] . multiplication sign (green) [T2], open circle (magenta) [T51 . 
filled diamond (blue) [14], filled square (red) [15], crossed dia¬ 
mond (cyan) [16], crossed square (blue) [T7: and open square 
(green) [18]. The SLAC data from the 1990’s are filled star 
(blue) [T9] and open diamond (magenta) [20]. The JLab data 
are asterisk (green) [21.! and filled triangle (blue) [22] . 



Fig. 4. Gm p /H p Gd extracted from cross section measurements 
versus Q 2 . The symbols are the same as in Fig. [3] Additional 
data points at the highest Q 2 , open square (magenta) [53] and 
open star (green) [54], were extracted from cross sections as¬ 
suming UpGep/Gmp = 1. The solid (dashed) line is a fit by 
Ref. [55] (Ref. [56]) described in the text. 


fits also included 31 normalization parameters for possi¬ 
ble systematic effects with cross sections measured with 
the different spectrometers and in different run periods. 
The fits which had a reduced y 2 < 1.16 had a maximum 
difference in their cross sections of 0.7%. 

At JLab, Ge p and Gm p were measured using the Rosen¬ 
bluth separation method at Q 2 between 0.4 to 5.5 GeV 2 as 
part of an experiment to measure inelastic cross sections 
on a range of nuclei m ■ The results are plotted in Figs. [3] 
and [4] as asterisk and agree with previous measurements. 

Instead of detecting the elastically scattered electron, 
(e,e'), an experiment which detected the elastically scat¬ 
tered proton to identify elastic reactions, (e,p), was run 
at JLab in 2002 [22] . The same experimental approach 
of extracting the form factors by measuring elastic cross 
sections at fixed Q 2 and different e by varying the beam 
energy was used. The experimental method takes advan¬ 
tage of the fact that the proton momentum is constant 
for all e at a fixed Q 2 . In addition for (e,p), the detected 
proton rate and the radiation corrections have a smaller 
dependence on e compared to (e, e') experiments. All this 
combines to reduce the e dependent systematic error com¬ 
pared to (e, e') experiments. The form factors were mea¬ 
sured at Q 2 = 2.64, 3.10 and 4.60 GeV 2 . In Fig. [3] and SI 
the measurements of Ge p and Gm p from the (e,p) reac¬ 
tion are plotted. The agreement between the form factors 
extracted by the different experiments is excellent. Ex¬ 
periments which detect either the scattered electrons or 
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the scattered protons have different systematics, so the 
agreement between the two techniques indicates that the 
experimental systematic errors are understood. 

With the success of the first JLab (e,p) experiment, a 
subsequent experiment, E05-017, was run at JLab in Hall 
C in 2007 [BT]. The experiment measured cross-sections at 
a total of 102 kinematic settings covering a wide Q 2 range 
from 0.4 to 5.76 GeV 2 with at least three e points per Q 2 . 
The emphasis was to measure at each Q 2 as wide an e 
range as possible. Fig. [5] plots the Q 2 versus e for all kine¬ 
matic points of E05-017. To obtain multiple e at the each 
Q 2 , 17 different beam energies were needed for the exper¬ 
iment. This number of beam energies in a relatively short 
time period demonstrate the amazing capabilities in the 
operation of the Continuous Electron Beam Accelerator at 
JLab. At Q 2 = 1 GeV 2 , thirteen e points were measured 
ranging from e = 0.05 to 0.98, with eight of the points 
above e = 0.8. Similarly for Q 2 = 2.3 GeV 2 , ten e points 
were measured ranging from e = 0.07 to 0.92, with five of 
the points above e = 0.7. The wide range of e at a fixed 
Q 2 allows a check of the non-linearity in the e dependence 
of the cross-section which would be a sign of two-photon 
exchange contributions effecting the cross-sections. The ef¬ 
fects from two-photon exchange contributions could have 
a dramatic e dependence near e = 1. 



Fig. 5. Kinematic coverage of the JLab experiment E05-017. 
Each kinematic point of the experiment is plotted for Q 2 ver¬ 
sus e. The solid and dashed lines are constant energy and the 
dotted lines are constant Q 2 . 


3.1.2 Neutron Form Factors 

The neutron has zero charge and therefore Geu has to 
be zero at Q 2 = 0. The slope of Ge u at Q 2 = 0 is 
related to the mean-square radius of the neutron (see 
Eq. (1381) for the proton). Also, the mean-square radius 
of the neutron, < r en 2 >, can be expressed in terms 
of the neutron-electron scattering length. The neutron- 
electron scattering length, b ne , can be determined from 


total transmission cross sections for epithermal and ther¬ 
mal neutrons scattering on the atomic electrons in no¬ 
ble gas targets [B2] and lead and bismuth targets mm 
E5]. Using results from the total transmission experiments, 
the Particle Data Group published a recommended value 
of < r en 2 >= -0.1161 ± 0.0022 fm 2 or dG En /dQ 2 = 
0.01935 ± 0.00037 fm 2 [55]. 

With no free neutron target, cross section experiments 
have to make measurements on a deuteron target to ex¬ 
tract the neutron form factors. Experiments have mea¬ 
sured cross sections for quasi-elastic single arm d(e, e') 
and coincidence d(e, e’p)n reactions. Primarily, the quasi¬ 
elastic ed reaction is a measurement of Gmu with lim¬ 
ited sensitivity to Gb„, since Geu is near zero and much 
smaller than the proton contribution to the cross section. 
Single arm quasi-elastic ed scattering by Ref. rrmrBniHH] 
can be used to extract Gmu but this requires theoreti¬ 
cal knowledge of the large final state interactions at low 
Q 2 , which leads to a sizable theoretical uncertainty. Co¬ 
incidence cross section measurements in the d(e, e'p)n re¬ 
action were done by Earn Detection of the neutron in 
coincidence reduces the theoretical uncertainty and the 
proton contribution in the extraction, but the uncertainty 
on the knowledge of the neutron detection efficiency be¬ 
comes important. In the 1990’s, Gmu was measured at 
MIT-Bates at Q 2 of 0.11, 0.18 and 0.26 GeV 2 using the 
d(e, e’n)p reaction and the data are plotted in Fig. E] as 
open squares. This experiment measured the neutron de¬ 
tection efficiency using the 2 H("/,pn) reaction. 

Extracting Gmu from the ratio of cross sections of the 
quasi-elastic d(e, e'n)p to d(e, e'p)n reactions is the least 
sensitive method to uncertainties in the calculation of the 
deuteron wave function, final state interactions and me¬ 
son exchange contributions. In the early 1970’s, at DESY, 
pioneering experiments measuring the ratio of quasi-free 
cross sections for the d(e, e'n) to d(e, e'p) reactions were 
performed by Ref. m and m- These experiments ex¬ 
tracted Gmu at Q 2 = 0.4, 0.57, 0.78, 1.0 and 1.5 GeV 2 
and the Gmu data are plotted in Fig. El as squares with 
cross. 

In 1995, at Bonn, the ratio of quasi-elastic d(e, e'n)p to 
d(e,e'p)n cross sections was used to extract Gmu- Gmu 
was measured at Q 2 of 0.13, 0.25, 0.42 and 0.61 GeV 2 
m and is plotted in Fig. El with empty diamond. The 
neutron detection efficiency was measured in situ using 
the 1 H('y,n~)n reaction. The photons were produced by 
bremsstrahlung in the hydrogen target. 

A series of measurements of the ratio of d(e, e'n)p to 
d(e, e'p)n cross sections were made in which the neutron 
detector efficiency was measured by taking the neutron de¬ 
tector to Paul Scherrer Institute (PSI) and using the neu¬ 
tron beam line. A LED system was used to monitor gain 
and baseline shifts at PSI and during the experiments. 
The first experiment was done at NIKHEF and Gmu was 
measured at Q 2 of 0.61 and 0.70 GeV 2 [75]. The next ex¬ 
periment was done at MAMI and measurements of Gmu 
were done at Q 2 of 0.24, 0.50, 0.65 and 0.78 GeV 2 [76]• The 
results from these two experiments are plotted as an aster¬ 
isk in Fig. El These were followed by more measurements 
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Fig. 6. Gmti/UtiGd versus Q 2 . The symbols corresponds to 
the data from open circle (green) [71], filled circle (black) [72], 
open square (black) |73 |. open diamond (black) [74], asterisk 
(blue) [75] and m, filled diamond (green) [77], open triangle 
(magenta) [?5j , square with cross inside (cyan) [IB] and dia¬ 
mond with cross inside (red) T9]. The solid line is a fit by 
Kelly [55] which was done in 2004 before the Hall B measure¬ 
ment [78] and excluded the data of Ref. m and m- 


of Gmu at MAMI which extend the Q 2 range. The latter 
experiment was done at Q 2 values of 0.071, 0.125, 0.36 
and 0.89 GeV 2 with statistical errors at the 1.5% level 
and are plotted in Fig. [ 6 ] as filled diamonds. Efficiency 
measurements at PSI were done before and after each ex¬ 
periment and the measurements were consistent. Within 
this series of experiments, the Gmu at matching Q 2 are in 
excellent agreement, but are smaller than Gmu measured 
at Bonn m and MIT-Bates m Ref. [ 8 T] suggested that 
the Bonn experiment had miscalculated their neutron effi¬ 
ciency, because a contribution from pion electroproduction 
was not taken into account when determining the neutron 
efficiency which led to an overestimate of Gmu- In their 
reply [S3) Bruins et al responded that in the peaking ap¬ 
proximation the contribution from pion electroproduction 
to their kinematics is negligible. In calculating the contri¬ 
bution from electroproduction, each paper uses a different 
data set to extrapolate to the measured kinematic region, 
so, as stated in Ref. [55], the only way to conclusively set¬ 
tle the disagreement is to measure pion electroproduction 
and photoproduction in the kinematics of the experiment. 

At the present time, the highest Q 2 measurement of 
G Mn was done at SLAC in the early 1980’s by measuring 
quasi-free ed cross sections m . This experiment measured 
Gmu from Q 2 of 2.5 to 10 GeV 2 . Another SLAC experi¬ 
ment [79 measured Gmu from Q 2 of 1.0 to 1.75 GeV 2 . In 


the 1990’s, at SLAC, a Rosenbluth separation experiment 
was done for the quasi-free ed reaction at Q 2 = 1.75, 2.5, 
3.25 and 4.0 GeV 2 and both Gmu and Geu were extracted 
E2- The Geu values were consistent with zero with large 
error bars. All experiments have consistent Gmu values in 
the region of overlapping kinematics and their Gmu values 
are plotted in Fig. [ 6 ] 

The experiment at JLab in Hall B using CLAS mea¬ 
sured Gmu in fine Q 2 bins from Q 2 between 1.0 to 4.8 GeV 2 
[78]. The data are plotted as open triangles in Fig. 

A unique feature of this experiment was a dual cell de¬ 
sign with liquid hydrogen and deuterium cells separated 
by 4.7 cm. This allowed measurement of the neutron de¬ 
tection efficiency by the H(e,e'n + )n reaction to be done 
simultaneously with cross section measurements. The Hall 
B data overlaps nicely with the SLAC measurements. Though 
at Q 2 near 4.8 GeV 2 , the Hall B data set suggest a less 
rapid Q 2 fall-off, then measured in the SLAC high Q 2 
data. 

Experiments have used elastic ed cross sections to de¬ 
termine the neutron form factors. The scattering by an 
electron from the spin 1 deuteron requires 3 form factors in 
the hadronic current operator, for the charge, quadrupole 
and magnetic distributions, Gc, Gq and GMd, respec¬ 
tively. In the original impulse approximation form of the 
cross section developed by Gourdin m, the elastic ed 
cross section is: 


da 

dfi 


% fA(Q 2 ) + H(Q 2 )tan 2 (|)), 
cLs* Mott \ Z J 


(26) 


where 

MQ~) = Gc(Q 2 ) + qV 2 Gq(Q 2 ) + 4^(1 +r])G 2 Md , 
B(Q 2 ) = ^rj(l + V ) 2 G 2 Md (Q 2 ), (27) 


with 77 = Q 2 /4M 2 . The charge, quadrupole and magnetic 
form factors can be written in terms of the isoscalar elec¬ 
tric and magnetic form factors as follows: 

° Md = JT p {GSmCs + \ G e C l ^ ( 28 ) 


where the coefficients Ce, Cq, Cl and Cs are Fourier 
transforms of specific combinations of the S- and D-state 
deuteron wave functions, u{r ) and w(r) [55]. The isoscalar 
(Gf) and isovector (Gf) magnetic (% = M) and electric 
(i = E) form factors are defined as : 


Gf — Gi p + Gi n and G\ — Gi P — G; n . (29) 


Theoretical knowledge of the deuteron wave function is 
needed to extract the form factors which is a major sys¬ 
tematic uncertainty. Both the magnetic and electric form 
factors can be extracted from the elastic ed cross sec¬ 
tions, though the Q 2 range is limited by the theoretical 
uncertainties. In the 1960’s, elastic ed cross section mea¬ 
surements were done which determined Gmu and G En for 
Q 2 < 0.12 GeV 2 [541155] . 
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Fig. 7. Geu extracted from the ed elastic reaction. Geti using 
the fit form of Eq. 11321) from Ref. [88] fitted to A(Q 2 ) data 
with either the Nijmegen solid line (black) or a Reid soft core 
dashed (red) NN potential in the theoretical calculation of the 
deuteron wave function. The dotted line (blue) is Eq. (1311) . The 
diamond points are from Ref. [89] . 


The 1971 DESY experiment of Galster et al. [SB] mea¬ 
sured elastic ed cross sections for forward scattered elec¬ 
trons for Q 2 up to 0.6 GeV 2 . At these kinematics, the 
cross section in Eq. (l26l) is dominated by the A(Q 2 ) term 
and GmcL contributes less than 5% to the A(Q 2 ) term. 
The A(Q 2 ) data was fitted using different deuteron wave 
functions and by using 

G Ep = ^A = ^ = Gn, (30) 

ftp H'n 

with different parametrization of G En . The lowest y 2 for a 
fit was obtained using the Feslibach-Lomon EZj deuteron 
wave function and the following fitting function: 

G En (Q 2 ) = - (y, G d (Q 2 ). (31) 

1 + 5.6t 

This fit is plotted in Fig. [7] as a dotted line. 

The most recent experiment to measure the elastic ed 
cross section to determine G En is that of Platchkov et al. 
[55]. These data extend to Q 2 of 0.7 GeV 2 , with signif¬ 
icantly smaller statistical uncertainties than all previous 
experiments. The form factor A(Q 2 ) is very sensitive to 
the deuteron wave function, and therefore to the NN in¬ 
teraction. Furthermore, the shape of A(Q 2 ) cannot be ex¬ 
plained by the impulse approximation alone. Corrections 
for meson exchange currents (MEC) and a small contri¬ 
bution from relativistic effects were found to significantly 
improve the agreement between calculations and the mea¬ 
sured shape of A(Q 2 ). When fitting the A(Q 2 ) data, a 
modified form of the Galster fit, 

G E „(Q 2 ) = -^^, (32) 


was used. Several NN potentials which including meson 
exchange currents as well as relativistic corrections were 
used to calculate the deuteron wave function. In Fig. [7] the 
fits of G En extracted from fitting a and b in Eq. m to 
the measured A(Q 2 ) are plotted when using the Nijmegen 
(black solid line) or a Reid soft core (red dashed) NN po¬ 
tential to calculate the deuteron wave function. Both fits 
to A(Q 2 ) had similar y 2 and the spread between the line 
gives a sense of the theoretical uncertainty in extracting 
G En from the elastic ed cross section. In 2001, an extrac¬ 
tion of G En was performed using the entire elastic ed cross 
section and polarization data [55] and the results are plot¬ 
ted in Fig. 0 as diamonds with the error bars showing the 
theoretical uncertainty. These data show the limit of using 
the ed elastic reaction to determine G En and the need to 
use the quasi-elastic ed polarization observables to extract 
Geu/Gmu which will be discussed in Sec. 13.2.21 

3.2 Double Polarization Experiments 

Both the recoil polarization method, and the asymmetry 
measurement using polarized target, have been used to 
measure the proton and the neutron form factors. Here we 
first describe the proton form factor results; the neutron 
form factor results will be discussed in the next subsection. 

3.2.1 Proton Form Factors 

The earliest polarization experiments, measuring the po¬ 
larization of the recoil proton [34] , or measuring the asym¬ 
metry using a polarized proton target [35 with unpolar¬ 
ized electron beams, were done to search for two photon 
effects. 

The first experiment with polarized electron beam and 
polarized target was done at the Stanford Linear Acceler¬ 
ator Center (SLAC) in 1970’s [50]. This experiment mea¬ 
sured the beam-target asymmetry A = at Q 2 = 

0.765 GeV 2 . The experiment showed that the results and 
the theoretical values were in good agreement if the signs 
of G Ep and Gmp are the same. 

The recoil polarization method was used for the first 
time in an experiment at the MIT-Bates laboratory to 
measure the proton form factor ratio Ge p /Gmp- This ex¬ 
periment determined G Ep /Gmp for a free proton [571138) . 
as well as for a bound proton in a deuterium target [39], 
at Q 2 -values of 0.38 and 0.5 GeV 2 . The success of this ex¬ 
periment highlighted the fact that the recoil polarization 
transfer technique would be of great interest for future 
measurements of G E and Gm at higher Q 2 values, for 
both the proton and the neutron. 

Next, using the same method of measuring the recoil 
polarization in 1 H(e, e'p) reaction, the ratio Ge p /Gm p was 
measured at MAMI at Q 2 -values of 0.373, 0.401 and 0.441 
Gev 2 [5TJ . The ratio results were found to be in agreement 
with those of Milbratlr et al. E3ED2 as well as Rosenbluth 
measurements. 

In the late 1990’s and 2000’s measurements using the 
recoil polarization method were made at JLab in Hall A 
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Fig. 8. Ratio h p Ge p /Gm p for Q 2 smaller than 1.7 GeV 2 
from number of different experiments from several laborato¬ 
ries; Open square (green) [371138] . filled diamond (black) [91] . 
filled circle (blue) [3j8j , filled square [92], filled triangle (black) 
[93] , crossed diamond (cyan) [94], open diamond (black) [951 . 
crossed square (magenta) [SB], open triangle (blue) [571, filled 
star (cyan) [98], open star (cyan) [99], open circle (red) [100 :. 
The curve is a 7 parameter fit given in Eq. 1441 with ratio 
p, p Ge p /Gm p constrained to 1 at Q 2 =0; the fit is of the Kelly 
type m, polynomial over polynomial, with 1/Q 2 behavior at 
large Q 2 . 


and Hall C [Mlf93l[Mll95ll96] at low Q 2 values, as cali¬ 
bration measurements for other polarization experiments. 
Two new high precision ratio measurements at low Q 2 
were made in Hall A at JLab; the first in 2006 measured 
the ratio in the range of Q 2 from 0.2 to 0.5 GeV 2 [57], the 
second in 2008 measured it at Q 2 of 0.3 to 0.7 GeV 2 [58] . 

The proton form factor ratio has also been obtained by 
measuring the beam-target asymmetry in the 1 H(e, e'p) 
reaction at a Q 2 of 1.51 GeV 2 in a Hall C experiment at 
JLab in elastic ep scattering [55]. This is the highest Q 2 
at which the Ge p /Gmp ratio has been obtained from a 
beam-target asymmetry measurement. The same method 
was used by the BLAST group at MIT-Bates [lOOj : this 
experiment measured the ratio at Q 2 values of 0.2 to 0.6 
GeV 2 with high precision. 

Figure [5] shows all the low Q 2 data obtained from re¬ 
coil polarization experiments [771IT1I81[5T1[5?1[581I5IH551[551 
HHHj and beam-target asymmetry measurements [5511 lOOj 
obtained at MIT-Bates, MAMI, and JLab. As can be seen 
from figure [51 data from different experiments are in gen¬ 
eral agreement. The slow decrease of the data starts at 
Q 2 « 0.5 GeV 2 and continues to 1.7 GeV 2 . 


A real break-through was made towards the under¬ 
standing of the internal structure of the proton, when two 
JLab Hall A and one Hall C experiments obtained the 
elastic electromagnetic form factor ratio of the proton, 
G p e /G p m at Q 2 ’s larger than 1 GeV 2 , from the measured 
recoil proton polarization components Pt and Pi , using 
the recoil polarization method. The first of these experi¬ 
ments measured the proton form factor ratios for Q 2 from 
0.5 to 3.5 GeV 2 in 1998 [7J[5], the second from 4.0, 4.8 and 
5.6 GeV 2 in 2000 [51110] and the third in 2007-8 up to 8.4 
GeV 2 mm . 

In the first JLab experiment GEp(l), elastic ep events 
were selected by detecting the scattered electrons and the 
recoiling protons in coincidence, using the two identical 
high-resolution spectrometers (HRS) of Hall A 11021 . One 
of the HRS was equipped with a focal plane polarime- 
ter (FPP) to detect the polarization of the recoil protons. 
The FPP consisted of two front detectors to track incident 
protons, followed by a graphite analyzer and two rear de¬ 
tectors to track scattered particles. The polarization of the 
recoiling proton was obtained from the asymmetry of the 
azimuthal distribution of the proton after re-scattering in 
the graphite analyzer of the polarimeter. 

In the second JLab experiment, GEp(2), the ratio, 
Ge p /Gm p was measured at Q 2 = 4.0, 4.8 and 5.6 GeV 2 
with an overlap point at Q 2 = 3.5 GeV 2 [51fTU[. Sev¬ 
eral changes were made compared to the first experiment, 
to extend the measurement to higher Q 2 . First, to in¬ 
crease the coefficient-of-merit (COM) of the focal plane 
polarimeter (FPP), a CH 2 analyzer was used instead of 
the graphite; hydrogen has much higher analyzing power 
[10311104] than carbon 11551 ; and, to increase the fraction 
of events with the proton interacting in the analyzer, the 
thickness of the analyzer was increased from 50 cm of 
graphite to 100 cm of CH 2 . Second, to achieve complete 
solid angle matching with the HRS detecting the proton 
and determining its polarization, a large frontal area lead- 
glass calorimeter was constructed and replaced the second 
HRS used in GEp(l). At the largest Q 2 of GEp(2) of 5.6 
GeV 2 , the solid angle of the electromagnetic calorimeter 
was 6 times that of the HRS. 

The results from the first two JLab experiments [71151 
151(15] . are plotted in Fig.[5]as the ratio p p Ge p /Gm p versus 
Q 2 , where they are compared with Rosenbluth separation 
data mmmmmmmm ■ As can be seen from this 
figure, for the polarization data at the larger Q 2, s the sta¬ 
tistical uncertainties are small, unlike those of the cross 
sections data, underlining the difficulties in obtaining Ge p 
by the Rosenbluth separation method at larger Q 2 ' s; the 
Rosenbluth data also show a large scatter among the re¬ 
sults from different experiments. 

The p p Ge p /Gm p ratio results from JLab [7][5][5] showed 
conclusively for the first time, a clear deviation of this ra¬ 
tio from unity, starting at Q 2 ~ 1 GeV 2 ; older data from 
mmm showed such a decreasing ratio, but with much 
larger statistical and systematic uncertainties, as seen in 
Fig. [5] The most important feature of the JLab data is 
the sharp decrease of the ratio p p Ge p /Gmp from 1, start¬ 
ing at Q 2 « 1 GeV 2 to a value of ~ 0.35 at Q 2 = 5.6 
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Fig. 9. The ratio fi p Ge p /Gm p from the two JLab experiments 
filled circle (blue) [7j[S] and filled square (red) [5] , compared to 
Rosenbluth separation results in (green), open diamond 1201 . 
open circle m, filled diamond [22], and open square EH 
I15II16II19| . The fit to the data is as in Gayou et al. [9]. 


GeV 2 , indicating that Ge p falls faster with increasing Q 2 
than Gm p , thus clearly highlighting a definite difference 
between the spatial distributions of charge and magneti¬ 
zation at short distances. This was the first definite ex¬ 
perimental indication that the Q 2 dependence of Ge p and 
Gm p are different. These results were very surprising at 
the time (1998-2002), as they appeared to contradict the 
previously accepted belief that the ratio h p Ge p /Gm p re¬ 
mains close to 1, a consensus based on the Rosenbluth 
separation results up to 6 GeV 2 , as illustrated in Fig. [9] 
As discussed above, the two methods available to de¬ 
termine the proton form factors Ge p and Gm p , the Rosen¬ 
bluth separation and polarization transfer, give definitively 
different results; the difference cannot be bridged by either 
simple re-normalization of the Rosenbluth data m, or by 
variation of the polarization data within the quoted sta¬ 
tistical and systematic uncertainties. This discrepancy has 
been known for sometime now, and has been the subject 
of extensive discussion and investigation. A possible ex¬ 
planation is the contribution from the hard two-photon 
exchange process, which affects the polarization transfer 
components at the level of only a few percent, but has 
drastic effects on the Rosenbluth separation results. This 
will be discussed in detail in section ITiil 

Following the unexpected results from the two first po¬ 
larization transfer experiments in Hall A at JLab, GEp(l) 
and GEp(2), a third experiment in Hall C, GEp(3), was 
carried out to extend the Q 2 -range to ~ 9 GeV 2 . Two new 
detectors were built to carry out this experiment: a large 



Q 2 (GeV 2 ) 


Fig. 10. All data for the ratio h p Ge p /Gm p obtained from the 
three large Q 2 recoil polarization experiments at JLab (filled 
circle (blue) HI) filled star (magenta) 11 i'H'i . filled square (red) 
[10] and filled triangle (black) 1101] '! compared to Rosenbluth 
separation data (green), open diamond |20| . open circle [2D . 
filled diamond HU- The curve is the same as in Figure [HJ a 7 
parameter fit given in Eq. 1441 


solid-angle electromagnetic calorimeter and a double focal 
plane polarimeter (FPP). The recoil protons were detected 
in the high momentum spectrometer (HMS) equipped with 
two new FPPs in series. The scattered electrons were de¬ 
tected in a new lead glass calorimeter (BigCal) built for 
this purpose out of 1744 glass bars, 4x4 cm 2 each, and a 
length of 20Ao, with a total frontal area of 2.6 m 2 which 
provided complete kinematical matching to the HMS solid 
angle. This experiment was completed in the spring of 
2008 and measured the form factor ratio at Q 2 of 5.2, 6.7 
and 8.5 GeV 2 . 

Figure [TO] shows the results from the three JLab ex¬ 
periments miwioiiMTon . as the ratio h p Ge p /Gm p versus 
Q 2 . The uncertainties shown for the recoil polarization 
data are statistical only. 

The striking feature of the results of the GEp(3) exper¬ 
iment is the continued, strong and almost linear decrease 
of the ratio with increasing Q 2 , albeit with some indica¬ 
tion of a slowdown at the highest Q 2 . The GEp(3) overlap 
point at 5.2 GeV 2 is in good agreement with the two sur¬ 
rounding points from the GEp(2) data [5MTD]. The GEp(3) 
experiment used a completely different apparatus in a Q 2 
range where direct comparison with the Hall A recoil po¬ 
larization results from the GEp(2) experiment is possi¬ 
ble. This comparison provides an important confirmation 
of the reproducibility of the results obtained with the re¬ 
coil polarization technique. Additionally, the results of the 
high-statistics survey of the e-dependence of Ge p /Gm p at 
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Q 2 = 2.5 GeV 2 , obtained from the GEp( 27 ) experiment 
m, which ran at the same time as the GEp(3) exper¬ 
iment is shown as a magenta star in Fig. 1101 and is in 
excellent agreement with the results from the GEp(l) ex¬ 
periment in Hall A [7J[§] at Q 2 = 2.47 GeV 2 . 

The results of the three JLab GEp experiments are 
the most precise measurements to date of the proton form 
factor ratio in this range of Q 2 , hence they represent a very 
significant advancement of the experimental knowledge of 
the structure of the nucleon. The proton electromagnetic 
form factor results from Jefferson Lab at high values of 
the four-momentum transfer Q 2 have had a big impact on 
progress in hadronic physics; these results have required 
a significant rethinking of nucleon structure which will be 
discussed in the theory section. 

3.2.2 Neutron Form Factors 

The early measurements of the form factors of the neutron 
are discussed in section 15.1.21 in this section only double 
polarization measurements are discussed. The recoil polar¬ 
ization and beam-target asymmetry, both techniques that 
have been used to measure Ge p and Gm p , also have been 
used to measure Geu and Gmu ■ However, as there are no 
free neutron targets, measurements of Geu and Gmu are 
more difficult than Ge p and Gm p - To make these mea¬ 
surements, complex light targets like 2 H and 3 He must 
be used in quasi elastic scattering. First, the recoil polar¬ 
ization experiments, and next the beam-target asymmetry 
experiments to extract Geu, will be described. 

The use of the recoil polarization technique to mea¬ 
sure the neutron charge form factor was made first at the 
MIT-Bates laboratory in the late 80’s using the exclu¬ 
sive 2 H(e, e'n)p reaction [ 107] . The advantage of using a 
deuteron target is that theoretical calculations predict the 
extracted neutron form factor results to be insensitive to 
effects like, final state interaction (FSI), meson exchange 
currents (MEC), isobar configurations (IC), and to the 
choice of the deuteron wave function |108U1091lll0l . In this 
experiment, the neutron form factor Geu was obtained 
from the measured transverse polarization component P/ 
of the recoiling neutron, and known beam polarization, P e , 
at a Q 2 of 0.255 GeV 2 . The relation between the polar¬ 
ization transfer coefficient Pt , the beam polarization, P e , 
and the measured neutron polarization component, P/, is 
P{ = P e Pt , the polarization transfer coefficient P t given 
by Eq. ©, is for a free neutron. This early experiment 
demonstrated the feasibility of extracting Geu from the 
quasi-elastic 2 H(e, e'n)p reaction with the recoil polariza¬ 
tion technique, with the possibility of extension to larger 
Q 2 values. 

The recoil polarization transfer method was next used 
at MAMI mum using the same reaction 2 H(e, e'n)p to 
determine Geu, at a Q 2 of 0.15 and 0.34 GeV 2 . How¬ 
ever, in this experiment the recoil neutron polarization 
components Pt and Pi were measured simultaneously, us¬ 
ing a dipole with vertical B-field to precess the neutron 
polarization in the reaction plane; the ratio P t /Pi, is re¬ 
lated directly to Ge/Gm as shown in Eq. (fill , again for a 


free neutron. As discussed earlier for the proton, the mea¬ 
surement of the ratio Pt/Pi, has some advantage over the 
measurement of Pt only; in the ratio the electron beam 
polarization and the polarimeter analyzing power cancel; 
as a result the systematic uncertainty is small. Also the 
model dependence for a bound neutron, which occurs via 
the dependence of the neutron wave function on the nu¬ 
clear binding, cancels in these polarization observables in 
leading order for the extraction of the form factor. 

Next, the electric form factor Geu was obtained at Q 2 
= 0.3, 0.6 and 0.8 GeV 2 from the measured ratio of po¬ 
larization transfer components, P t /Pi , in another exper¬ 
iment at MAMI EH; Glazier et al. concluded that the 
results from this experiment were in good agreement with 
all other Geu results from double-polarization measure¬ 
ments. 

The first double polarization experiment to measure 
GEn at JLab by Madey et al. |11211113j . obtained the neu¬ 
tron form factor ratios Ge/Gm at Q 2 values of 0.45, 1.13 
and 1.45 GeV 2 using the same method of measuring the 
recoil neutron polarization components Pt and Pi simulta¬ 
neously. The neutron charge form factor Geu was calcu¬ 
lated from the measured ratio, using the best-fit values of 
G Mn - This was the first experiment that determined Geu 
with small statistical and systematic uncertainty to rela¬ 
tively high Q 2 . Madey et al. concluded that a successful 
theoretical models must be able to predict both neutron 
and proton electromagnetic form factors simultaneously. 
The neutron electric form factor is more sensitive to small 
components of the nucleon wave function, and differences 
between model predictions for Geu tend to differ with in¬ 
creasing Q 2 ; hence the new data from this experiment to 
larger Q 2 provided a challenging test for theoretical model 
calculations. 

The first measurement of Geu using the beam-target 
asymmetry was made at NIKHEF at a Q 2 of 0.21 GeV 2 
using the 2 H(e, e'n)p reaction [40j . The experiment used a 
polarized electron beam from the storage ring and a vec¬ 
tor polarized deuterium gas target, internal to the stor¬ 
age ring; Geu was extracted from the measured sideways 
spin-correlation parameter in quasi-free scattering. Pass- 
chier et al. concluded that their result puts strong con¬ 
straints on Geu up to Q 2 = 0.7 GeV 2 when combined 
with the measured Geu slope from Kopecky et al. |124) 
at Q 2 =0 GeV 2 and the elastic electron-deuteron scatter¬ 
ing data from Platchkov et al. [55] . 

The neutron electric form factor at Q 2 = 0.5 and 1.0 
GeV 2 was extracted from measurements of the beam-target 
asymmetry using the 2 H(e, e'n)p reaction in quasi elastic 
kinematics, at JLab in Hall C [511152] : the polarized elec¬ 
trons were scattered off a solid polarized deuterated am¬ 
monia (ND 3 ) target in which the deuteron polarization 
was perpendicular to the momentum transfer. This was 
the first experiment to obtain Geu at a relatively large 
Q 2 using a polarized target. 

There was a measurement of Geu/Gmu at the MIT- 
Bates lab in mid 2000’s, using a longitudinally polarized 
electron beam and a vector-polarized 2 H target internal to 
the storage ring over a range of Q 2 between 0.10 and 0.55 
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Fig. 11. The complete data base for Geu/Gm n from double 
polarization experiments; recoil polarization with deuterium 
target |551I41U42U111U112U113| . asymmetry with polarized deu¬ 
terium target: and asymmetry with polarized 

bols for the data are filled square (red) |1071 . filled diamond 
(red) [42], open triangle (red) [411 . multiplication sign (red) 
m , crossed diamond (red) [lllj . open star (blue) HIT] , filled 
star (blue) m, open diamond (blue) 11211 . filled circle (ma¬ 
genta) [1201 . asterisk (black) |112| . plus sign (black) [32], 
crossed diamond (green) 11141 . filled triangle (red) [321, cross 
diamond (red) 1122 . The solid (blue) line is a fit using Eq. (1431) . 
Geti/Gm-h from the DSE model of Ref. |123 | are plotted as a 
dash-dotted (cyan) line. 


GeV 2 [114] : in this experiment the quasi-elastically scat¬ 
tered electrons were detected in coincidence with recoil 
neutrons in the BLAST detector. They used the parametriza- 
tion of Friedrich and Walcher 125} for Gmu to calculate 
Geu ■ Geis et al. | I114] concluded that their data are in ex¬ 
cellent agreement with VMD based models of Lornon m 
and Belushkin [127) , and also agree with the meson-cloud 
calculation of Miller D2S1- 

All the experiments described above used either a po¬ 
larized or an unpolarized deuterium target. 

In 1984 Blankleider and Woloshyn [129] . proposed that 
to measure Geu or Gmu , a polarized 3 He target could be 
used. Their argument was that the ground state of 3 He 
is dominated by the spatially symmetric S-state in which 
the two proton spins point in opposite directions, hence 
the spin of the nucleus is largely carried by the neutron. 
Therefore, polarized 3 He target effectively serve as a po¬ 
larized neutron target; and in the quasi-elastic scattering 
region the spin-dependent properties are dominated by the 
neutron in the 3 He target. 

The first two experiments that used a polarized 3 He 
target and measured the asymmetry with polarized elec¬ 


trons in spin-dependent quasi-elastic scattering were done 
at MIT-Bates Laboratory [11511116] ; these experiments ex¬ 
tracted the value of Geu at a Q 2 =0.16 and 0.2 GeV 2 , using 
the model of Blankleider and Woloshyn m- However, 
Thompson et al. [[ 1161 pointed out that significant correc¬ 
tions were necessary at Q 2 = 0.2 GeV 2 , for spin-dependent 
quasi elastic scattering on polarized 3 He according to the 
calculation of Laget IllOj ; hence no useful information on 
GEn could be extracted from these measurements; but 
Thompson et al. m concluded that at higher Q 2 values 
the relative contribution of the polarized protons becomes 
significantly less and a precise measurements of GEn using 
polarized 3 He target will become possible. 

The neutron electric form factor Geu was obtained in 
the early 1990’s in several experiments at MAMI; these ex¬ 
periments measured the beam-target asymmetry in the ex¬ 
clusive quasi-elastic scattering of polarized electrons from 
polarized 3 He in the 3 He(e, e'n)pp reaction [11711118111191 
EM- The first of these experiments at MAMI obtained 
Geu at Q 2 = 0.31 GeV 2 EZI; and the next experiment 
measured Geu at Q 2 of 0.35 GeV 2 [ 118] and 0.67 GeV 2 
[11911120] using the same experimental setup. The value 
of Geu at Q 2 of 0.35 GeV 2 m was later corrected by 
Golak et al. m, based on Faddeev solutions and with 
some MEC corrections. The size of these corrections is 
expected to decrease with Q 2 , although the corrections 
become increasingly difficult to calculate with increasing 

Q 2 . 

The GEn(l) experiment in Hall A at JLab measured 
the ratio Geu/Gmu in 2006 at a Q 2 = 1.72, 2.48, and 
3.41 GeV 2 using the reaction 3 He(e, e'n)pp in quasi-elastic 
kinematics [33] . Longitudinally polarized electrons were 
scattered off a polarized target in which the nuclear po¬ 
larization was oriented perpendicular to the momentum 
transfer. The scattered electrons were detected in a mag¬ 
netic spectrometer in coincidence with knocked out neu¬ 
trons, that were detected in a large hadron detector. The 
ratio Geu/Gmu was obtained from the measured beam- 
target asymmetry. Riordan et al. [33] concluded that this 
-experiment more than doubled the Q 2 range over which 
Geu is known, and this fact greatly sharpens the mapping 
of the nucleon’s constituents, and provides a new bench¬ 
mark for comparison with theory. 

There is yet another recent experiment at MAMI, that 
measured the charge form factor of the neutron at a Q 2 of 
1.58 GeV 2 using the polarized 3 He target and longitudi¬ 
nally polarized electron beam m ■ To reduce systematic 
errors, data were taken for four different target polariza¬ 
tion orientations. The data of this experiment are in very 
good agreement with the data of Riordan et al. [33]. 

Figure fill shows the ratio Geu/Gmu versus Q 2 ob¬ 
tained from all double polarization experiments in the last 
two decades. The results from each experiment are shown 
as different symbols as explained in the figure caption. 
The solid curve is a fit to the data as given by Eq. m- 

The first experiment to obtain the magnetic form fac¬ 
tor of the neutron, Gmu , from polarization observables 
was done at the MIT-Bates laboratory. This experiment 
obtained Gmu from the measured beam-target asymme- 
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Fig. 12. Gmti/UtiGd versus Q 2 . The data from polarization 
experiments: the open star (magenta) 130| . and the filled 
square (cyan) |1331 . The data from cross sections (black) are 
same as in Figure [S] The solid line is a fit by Kelly |55) . 


try in inclusive quasi-elastic scattering of polarized elec¬ 
trons from polarized 3 He target at Q 2 of 0.19 GeV 2 |130j : 
the uncertainty on Gmu was dominated by the statistics, 
with a relatively small contribution from model depen¬ 
dence of the analysis. The second experiment was done 
at JLab in Hall A; this experiment extracted Gmu for 
Q 2 values between 0.1 and 0.6 GeV 2 , by measuring the 


transverse asymmetry in the 3 He(e, e') reaction in quasi- 
free kinematics [13111132111331 . The values of Gmu were 
obtained with a full Faddeev calculation at Q 2 of 0.1 and 
0.2 GeV 2 , and in the plane wave impulse approximation 
(PWIA) at Q 2 of 0.3 to 0.6 GeV 2 . It was argued by the 
authors of this paper that the PWIA extraction of Gmu 
is reasonably reliable in the Q 2 range of 0.3 to 0.6 GeV 2 ; 
however, a more precise extraction of Gmu requires fully 
relativistic three-body calculations. The Gmu data from 
both double polarization experiments |130lll31lll32lll33| 
are shown in Fig. [12l together with the data from ear- 

There is some scatter in the data, but they are 
close to Gmu/HuGd ~ 1.0, except two open circles data 
point at a larger Q 2 values. 


3.3 Discrepancies between Cross Section and 
Polarization Results 

As discussed in Sec. 13.2. H and illustrated in Figs. l9l and flOl 
there is a clear discrepancy between the extracted values 
of the proton form factor ratio, Ge p /Gmp, from polar¬ 
ization transfer and double polarization experiments, and 


those obtained from cross section measurements using the 
Rosenbluth separation technique. As also mentioned pre¬ 
viously, one possible explanation for this discrepancy is 
related to the hard two-photon exchange process; more 
properly, it is the interference between single-photon and 
two-photon exchange processes that has the potential to 
complicate the extraction of the form factors in Rosen¬ 
bluth separation experiments. 

Two groups simultaneously suggested that the differ¬ 
ence between cross section and double polarization results 
might be due to previously neglected two-liard-photon ex¬ 
change processes; these were Guichon and Vanderhaeghen 
[50], and Blunden et al. m- In general, cross section 
data require large radiative corrections, whereas double¬ 
polarization ratios do not. This is because radiative cor¬ 
rections affect the longitudinal and transverse polarization 
observables similarly, and thus the residual correction for 
double polarization is at the few percent level. A number 
of calculations of the two-hard-photon contribution have 
been published over the last decade. A partial list of calcu¬ 
lations of the contribution of the two-hard-photon process 
to the cross section includes ICTlll35Ul^Slll37l[T38H139l . 

The signature of two-photon exchange processes would 
be an observed e-dependence; in the case of cross section 
experiments the effect is an non-linearity in the reduced 
cross section as a function of e, whereas in the polarization 
transfer experiments, one might expect an e-dependence in 
the ratio of the measured proton polarization components. 
There have been several experimental attempts over the 
last decade to search for these dependencies. 


3.3.1 Results from the Ge p { 2y) Experiment 

The Ge p {^"i) experiment in Hall C at JLab recently mea¬ 
sured the ratio — \/(l + e)/2e(P t /Pi), which strictly equals 
Ge p /Gm p in the Born approximation, at a central value 
Q 2 =2.49 GeV 2 , and for three values of e: 0.152, 0.635 
and 0.785, with very small error bars [106]. P t and Pt 
are the transverse and longitudinal components of the po¬ 
larization transferred to the proton, respectively. Simul¬ 
taneously, values of Pi/Pi,Bom were obtained at the two 
larger e values, using the lowest e data point to determine 
the analyzing power of the polarinreter for the common 
central proton momentum of 2.06 GeV/c. For these data, 
radiative corrections were calculated using the model in¬ 
dependent calculation of P323- The calculation includes 
the vacuum polarization correction, the electron vertex 
correction, and internal and external bremsstrahlung cor¬ 
rections. In general, polarization observables are rather 
insensitive to radiative corrections, as they are by defini¬ 
tion a ratio of a polarized cross section to an unpolarized 
one. Moreover, the bremsstrahlung correction can be re¬ 
duced drastically by applying missing mass or inelasticity 
cut, as was done in this experiment through elastic event 
selection. 

The results of these two m easuremen ts are shown in 
Fig. US] The measured ratio — \/(l + e)/2e(P t /Pi) appears 
to have no e dependence within the small statistical and 
systematic uncertainties of the experiment. In contrast, 
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Fig. 13. Results of the Ge p ( 2y) experiment at JLab i 106] . 
The ratio plotted in the top panel is equal to p v Ge p /Gm p in 
the Born (one-photon-exchange) approximation with the data 
of Ref. 11061 plotted as filled circles (blue) and the data point 
of Ref. [8] plotted as empty triangle (black). The theoretical 
curves are: long dot-dashed (magenta) 11411 . short dot-dashed 
(green) |142| . long dashed (cyan) 1 14311144| . solid (black) [ST] , 
dotted (red) 11451 with wave functions from m, and short- 
dashed (red) [145] with wave functions from Ml- The bottom 
panel shows Pe/Pe,Born for the two higher e points of the exper¬ 
iment plotted as filled circles (red). The lowest e point was used 
to determine the polarimeter analyzing power. The curves are 
fits from HE]; the dashed (red) curve corresponds to ”Fitl” 
and the solid (blue) curve corresponds to ”Fit2” as specified 
in this reference. 


the ratio Pi/Pi,Bom, displayed in the bottom panel of 
Fig. 1131 shows a systematic deviation from unity at the 
largest e value of up to 4.5 standard deviations. Such a 
behavior can be explained (see the curve in lower panel 
of Fig. n~3l) within the context of recent work described in 
(145] which shows that the corrections to the three form 
factors required in the presence of the interference of the 
one- and two-photon terms do not cancel one another as 
e -1 1 . 


3.3.2 Results from e + /e Scattering Experiments 


The most direct way to characterize a hard two-photon 
contribution to the elastic ep cross section is to compare 
e + p and e~p scattering. There are recent results from two 
experiments (VEPP-3 in Novosibirsk and in Hall B at 
JLab) which attempt to determine the two-hard-photon 
contribution via measurements of the ratio, R, of the elas¬ 
tic e + p and e~p scattering cross sections. In addition, the 
OLYMPUS experiment at DESY (which also will measure 
this ratio) is currently in the data analysis phase, and re¬ 
sults are expected to be published soon. 

In general, the lepton-proton elastic scattering cross 
section is proportional to the square of the sum of the 
Born amplitude and all higher-order QED correction am¬ 
plitudes. The ratio of e ± p elastic scattering cross sections 
can be written [148] as: 

_ Cr(e + p) ^14 Seven J 27 Sbrem 

(j(e p) 1 T Seven T J27 T Sbr em 
~ 1 2 (^ 2 ^ + fibrem )/ (1 “h even ) , (33) 


where 5 eV en is the total charge-even (relative to lowest or¬ 
der) radiative correction factor, and J 27 and Sbrem are the 
fractional two-photon-exchange and lepton-proton bremsstrahlung 
interference contributions, respectively. After calculating 
and correcting for the charge-odd Sbrem term, the cor¬ 
rected cross section ratio is: 


R' 


1 - 


25■ 


27 


(1 T Seven ) 


(34) 


The hard two-photon contribution to the ep scattering 
cross section is: 


J27 — 


2 Re(M\ 1 M h 2 “ rd ) 


| M 


17I 


(35) 


where AJi 7 and are the single- and two-photon 

hard scattering amplitudes, respectively. For the purposes 
of comparing to theoretical predictions, the results are 
sometimes presented as the ratio R 27 = (1 —$ 27 )/( 1+^27 )• 
The results for i? 2 7 from Runs I and II at VEPP-3 
[149] are shown in Fig. [ll] together with several theoreti¬ 
cal predictions. At JLab, in Hall B, the CLAS collabora¬ 
tion obtained two-photon exchange data for Q 2 between 
0.5 and 3 GeV 2 , with 0.15 < e < 0.95 (148j . Results from 
this experiment for the ratio R' are shown in Fig. [15] as 
a function of both Q 2 and e. The data from each of these 
experiments indicate that the hard two-photon-exchange 
effect is significant, and they are in moderate agreement 
with several two-photon-exchange predictions which also 
explain the form factor ratio discrepancy at higher Q 2 
values, thus pointing to two-photon-excliange as a likely 
source of at least part of the discrepancy. 


3.4 The Proton Charge Radius 

Non-relativistically, the elastic ep cross section is related 
to the product of the Mott cross section for a point-like 
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Fig. 14. Experimental data together with theoretical predic¬ 
tions for the ratio i? 2 7 as a function of t or Q 2 . The top and 
bottom panels correspond respectively to Run-I and Run-II. 
The data points are: open squares (green) loll: , downward tri¬ 
angles (green) [151 1. diamonds (green) 11521 . upward triangles 
(green) 153] . and circles (black) [1491 . The theoretical curves 
are from 11541 dash-dotted (green), 11421 thin solid (red), [601 
thick solid (blue), tl'55] long-dashed (black), 11561 short-dashed 
(magenta), and [157] dotted (black). 



Fig. 15. The ratio of e + p/e p cross sections corrected for 
Sbrem as a function of e at Q 2 = 1.45 GeV 2 (top) and as a 
function of Q 2 at e = 0.88 (bottom). The filled circles (blue) 
are from H18| . The inner error bars are the statistical uncer¬ 
tainties and the outer error bars are the statistical, system¬ 
atic and radiative-correction uncertainties added in quadra¬ 
ture. The line at R' = 1 is the limit of no two-photon-exchange 
contribution. The theory curves are: dotted line (black) - [142] . 
solid (magenta) and dashed (red) - 11581 including N only and 
N + A intermediate states, respectively, dot-dashed line (cyan) 
- [eh]. The open circles (green) show the previous world data 
(at Q 2 > 1 GeV 2 for the top plot) ' 135| . 


spin-i electron and the Fourier transform of the charge 
and/or magnetization density of the target nucleon, as 
follows: 


from which it follows that 


< r Ep 2 >= -6 


dG Ep 
dQ 2 


atQ 2 —0 


(39) 


@(de) — &Mott X 


p(r)e iqr d 3 


I volume 


(36) 


where p(r) is either the electric or the magnetic spatial 
distribution function. For the particular case of the electric 
form factor G Ep (Q 2 ), it follows that for short distances it 
can be expanded in terms of even moments of the distance 


< i% > as: 


G Ep = 1 — -Q 2 < r Ep 2 > + < r Ep 4 >••• 


(37) 


Hence, for very small distances within the nucleon, the 
mean-square radius of the proton can be obtained from 
the derivative of Eq. 


dG Ep 1 I 21 
~dQ 2" = _ 6 ' TEp la «2 2=0 


(38) 


Similar relations hold for the magnetic form factor, Gm p (Q 2 ), 
and the magnetic radius, < i' m p 2 >• The cumulative cross 
section data from electron scattering experiments at low 
Q 2 have been used to obtain values of < r Ep 2 > [160111611 
1162U16311164J . The extraction of < rjf p 2 > is more diffi¬ 
cult to obtain as its contribution to the cross section is 
suppressed by the factor r (see Eq. |0). 

In a completely complementary fashion, the proton ra¬ 
dius can also be obtained from precise measurements of 
the Lamb shift energies either in the hydrogen atom m 
or in muonic hydrogen [166] : indeed, a (different) fraction 
of the Lamb shift is related to the finite size of the pro¬ 
ton nucleus in each case. In both cases, the Lamb shift 
under consideration depends upon the overlap of the ap¬ 
propriate S-, P-, or D-state wave functions and the proton 
nucleus; it is in this way that the finite proton size con¬ 
tributes to the Lamb shift. Recent measurements of the 
muonic Lamb shift energies at PSI have produced values 



























































V. Punjabi, C.F. Perdrisat, M.K. Jones, E.J. Brash, and C.E. Carlson: The Structure of the Nucleon 


19 


of < te p 2 > which are smaller than the mean value of all 
electron scattering experiments (the so-called CODATA 
value - see Refs. [167U1681I169I 1 by about 4% (a 7 a dif¬ 
ference). This discrepancy has become known as the pro¬ 
ton radius puzzle, and its resolution has become a topic 
of great current interest, and the aim of several new and 
novel experimental efforts. We focus here on the experi¬ 
mental determinations of the proton radius; indeed, there 
has been a plethora of works concerned with the theoreti¬ 
cal aspects of the problem. We refer the reader to [6] and 
m, and references therein, for further details. 


3.4.1 Previous Results 

In Fig. [TH we show several determinations of the RMS 
proton charge radius, r p = ^< rE P 2 > over the last sev¬ 
eral decades. Earlier extractions were based on data from 
elastic ep scattering experiments at Orsay im, Stan¬ 
ford Ena, Saskatoon m and Mainz [17U174J . together 
with the various re-analyses of these world data II6OIII6II 
11621116311164] . In general, the extraction of the proton ra¬ 
dius from these data involves fitting the form factors with 
either simple mathematical or in some cases theoretically 
inspired parametrizations in order to determine the ra¬ 
dius. As can be seen clearly in Fig. [TGl the consistency 
between various approaches has improved over the years. 
However, several issues related to the extraction procedure 
remain. There is, inherently, a model dependence uncer¬ 
tainty in the extraction (as evidenced for example in the 
discrepancy between the open circles and open diamonds 
in Fig. na but this is typically not included in the quoted 
uncertainty for a given extraction. Moreover, it appears at 
this point that the treatment of systematic uncertainties 
in many experiments was overly optimistic. This can be 
inferred from the fact that in many global fits, the y 2 per 
degree of freedom is larger than unity. In some analyses, 
the absolute normalization of the cross section data was 
allowed to float, but this is done at the expense of sensi¬ 
tivity to the radius. In addition, there are issues connected 
with the range of Q 2 data that is included in the global 
fits. It can be seen from Eq. (EZD that the coefficients of 
successive < r)7 > terms increase with order. Therefore, 
it is not possible to define a value of Q 2 where any one 
term in the expansion sufficiently dominates such that its 
value could be fixed and then used in fitting to lower Q 2 
data. 

In order to address some of the issues related to sys¬ 
tematic uncertainties in the cross section data, the Mainz 
A1 collaboration at MAMI [59] measured 1422 precise rel¬ 
ative ep cross sections in the low Q' 2 region (0.0038 GeV 2 
to 0.98 GeV 2 ) and a wide range of beam energies and 
scattering angles. Experimental systematic uncertainties 
were controlled by using one spectrometer as a luminos¬ 
ity monitor, and then moving the other two spectrom¬ 
eters through multiple, overlapping angle settings. The 
cross section data were subsequently fit with a variety 
of functional forms, in order to assess model dependent 
uncertainties. Interestingly, it was found that satisfactory 
goodness of fit could only be obtained through the use of 
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Fig. 16. Proton radius determinations over time. Electronic 
measurements seem to settle around r p =0.88 fm, whereas the 
muonic hydrogen value is at 0.84 fm. The solid circles (blue) are 
all electron scattering measurements; Orsay 1962 nza. Stan¬ 
ford 1963 |172| . Saskatoon 1975 |173l . Borkowski 1975 [lT], and 
Mainz 1980 11741 are early measurements, whereas the more re¬ 
cent measurements are Mainz 2010 [59] and JLab 2011 I175 1 . 
The open circles (blue) Il60lll61lll62l and open diamonds 
(green) |1631I164] denote reanalyses of the world electron scat¬ 
tering data using various fitting functions (See Sec. 14.1.Hi . The 
solid triangles (red), (CODATA) originate from laser spec¬ 
troscopy of atomic hydrogen and advances in hydrogen QED 
theory Il67lll68lll69l . The solid squares (black) are the values 
extracted from muonic hydrogen [77511777] , 
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more flexible mathematical fitting functions (polynomials 
or splines) as opposed the more traditional physically mo¬ 
tivated forms, such as dipoles. The final extracted value for 
the proton radius was r p = 0.879 ± 0.005 stat ± 0.004 S!/St ± 
0.002 mo dei ± 0.004 ffrollp , where the final uncertainty comes 
from the polynomial vs. spline difference. 

In parallel with the cross section measurements, the 
Jefferson Lab LEDEX collaboration i 17311178ii f?9 mea¬ 
sured the proton form factor ratio p p Ge p /Gmp using po¬ 
larization transfer in the Q 2 ~ 0.1 —> 0.4 GeV 2 region. 
It is particularly interesting that while measurements of 
the form factor ratio do not give the proton radius di¬ 
rectly, accurate and precise knowledge of the ratio helps 
to constrain normalizations of cross section data during 
fits, which in turn leads to an improved value of the ex¬ 
tracted radius. The analysis of Ref. m gives r p = 0.875 
± 0.008 exp ± 0.006 fa, in agreement with the value ex¬ 
tracted in the Mainz analysis. It is interesting to note that 
this analysis uses an entirely independent data set from 
the Mainz analysis. 

The Lamb shifts in both normal electronic as well as 
muonic hydrogen are sensitive to the non-point-like struc¬ 
ture of the proton nucleus. While this effect is indeed mi¬ 
nuscule, the ability of quantum electrodynamics (QED) to 
predict the energy levels of hydrogen with remarkable ac¬ 
curacy, together with experimental measurements of the 
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Fig. 17. Proton charge radii r p obtained from hydrogen spec¬ 
troscopy. The proton charge radius can best be extracted from 
a combination of the 1S-2S transition frequency 1180) and one 
of the indicated 2S-NLj transitions nsntrssi . The thin ver¬ 
tical line (cyan) is the weighted average of the data points 
shown; the error bar on this average is represented by the ver¬ 
tical dashed lines. The value from muonic hydrogen TMT77] 
is represented by the thick vertical line (red); the thickness of 
this line represents the error bar. 


relevant transition frequencies to around one part in 10 11 , 
results in a measurable sensitivity to the proton radius. 
The energies of the S-states in hydrogen are given by: 

E(nS)~-^ + ^- (40) 

n z n a 

where n is the principal quantum number, and Lis de¬ 
notes the Lamb shift of the IS ground state which is given 
by QED and contains the effect of the proton charge ra¬ 
dius, r p . Numerically, Lis — (8172 + 1.56 r p (fm) 2 ) MHz, 
so the finite size effect on the IS level in hydrogen is about 
1.2 MHz. The different n-dependence of the two terms in 
Eq. (l40ll . permits the determination of both Roo and r p 
from at least two transition frequencies in hydrogen. The 
approach taken to date has been to use the 1S-2S tran¬ 
sition m, which has been measured to a stunning ac¬ 
curacy of one part in 10 15 , together with one of the 2S- 
NLj transitions [18111182] . The former transition is max¬ 
imally sensitive to r p , whereas the latter contain much 
smaller Lamb shift contributions, due to the 1/n 3 scaling 
in Eq. (l40l) . 

In Fig.[T7| we show the values of r p extracted from the 
various transition combinations. The r p values extracted 
from measurements in normal electronic hydrogen favor a 
r p value of ~ 0.88 fm, consistent with the world-averaged 
electron scattering results. The discrepancy between the 
combined value from just electronic hydrogen alone, as ob¬ 
tained in the elaborate CODATA adjustment of the fun¬ 
damental constants \m\ . and the muonic hydrogen value, 
is about 4.4cr. 

The difference between the muonic hydrogen deter¬ 
mination of the proton radius and the results obtained 
from either electron scattering or from electronic hydro¬ 
gen transitions is highly enigmatic. The theoretical QED 


calculations of the portion of the Lamb shift that is not 
due to the finite proton size have been checked and re¬ 
checked by multiple independent groups. Also, the many 
measurements of the transition frequencies in electronic 
hydrogen are in agreement with one another. Thus, from 
an experimental point of view, the solution to this puzzle 
is not obvious at this time. On the theoretical side, there 
has been a significant amount of interest in this problem 
as well; this is discussed in more detail in Section [I] 


3.4.2 Future Experiments 

Experimental approaches to unraveling the proton radius 
puzzle break down into several distinct categories: improv¬ 
ing the precision of the proton radius extraction from elec¬ 
tron scattering experiments, extending the atomic spec¬ 
troscopy measurements to other ions or exotic atoms, and/or 
determining the proton radius in elastic muon-proton scat¬ 
tering. 

At Jefferson Laboratory, in Hall B, there is an ap¬ 
proved experiment which aims to improve the electron 
scattering radius determination by extending the Q 2 range 
from the 0.0038 GeV 2 of the Mainz experiment down to 
1-2 x 10~ 4 GeV 2 [183] . This is indeed a challenging ex¬ 
periment. If the proton radius is to be determined in this 
experiment with similar accuracy compared to the cur¬ 
rent world-averaged result from electron scattering, then 
relative cross sections need to be determined at the 0.2% 
level. Fortunately, the large cross sections and correspond¬ 
ing event rates for low Q 2 scattering make statistical un¬ 
certainties of 0.1% achievable. However, systematic uncer¬ 
tainties must also be controlled at a similar level. One of 
the largest sources of systematic error is the determina¬ 
tion of the electron scattering angle; at the lowest Q 2 value 
of the experiment, the scattering angle is only 10 rnrad, 
and due to the severe angular dependence of the Mott 
cross-section at small angles, a 10 /zrad knowledge of the 
scattering angle is needed to limit shifts in the relative 
cross sections to 0.2%. 

At MAMI, an experiment is underway that aims to 
measure the proton electromagnetic form factors in ep 
scattering at very low momentum transfers by using a 
technique based on initial state radiation m ■ The basic 
premise is that initial state radiation degrades the energy 
of the incoming electron so that the momentum transfer 
to the proton can be quite low. The outgoing electron an¬ 
gle and energy are measured as usual, and together with 
theoretical input, an accurate form factor can in principle 
be obtained at Q 2 values as low as 10 -4 GeV 2 . The full 
experiment ran in 2013, and the analysis is continuing. 

As can be seen from Eq. (H3D, the determination of 
r p using transitions in hydrogen relies heavily on a pre¬ 
cise determination of Roo- Currently, Roc is known from 
measurements of transition frequencies in hydrogen and 
deuterium. However, since the correlation between r p and 
Roc is significant, new precise measurements of Roo (at 
the level of a few parts in 10 12 ) that are independent of 
r p could potentially impact the extracted value of r p from 
electronic hydrogen measurements. There are several such 




























V. Punjabi, C.F. Perdrisat, M.K. Jones, E.J. Brash, and C.E. Carlson: The Structure of the Nucleon 


21 


efforts currently underway. These include single- m and 
two-photon m measurements in hydrogen, as well as ex¬ 
periments which aim to use laser spectroscopy of neutral 
and ionic helium 


Noticeably, a crucial missing piece of the proton radius 
puzzle is a measurement of the proton radius using muon 
scattering; indeed, there is a proposal at PSI, the MUSE 
Experiment, to make just such a measurement |190] . The 
MUSE experiment will measure elastic pp scattering to a 
minimum Q 2 of 0.002 GeV 2 - about half the lower limit 
of the Mainz electron scattering experiment - using both 
positively and negative charged incident muons, so that 
any possible two-photon effects can be taken into account 
directly from the data, rather than relying on theoreti¬ 
cal calculations. In addition, the experiment will simulta¬ 
neously collect ep scattering data, so that the extracted 
proton radius from muon and electron scattering can be 
compared directly within a single experimental appara¬ 
tus. Preliminary estimates are that the proton radius that 
can be extracted from muon scattering will be similar in 
precision to that extracted from the Mainz experiment. 


3.5 Flavor Separation of Nucleon Form Factors 

Charge symmetry implies that the proton and neutron 
wave functions are identical under the interchange of the 
up and down quark contributions. Measurements of asym¬ 
metries in parity non-conserving electron scattering on the 
proton have found that the strange quark form factors are 
small (see the review article by Ref. [44]). Ignoring the con¬ 
tributions of higher mass quarks, the proton and neutron 
form factors can be written in terms of the contributions 
from the up and down dressed quark form factors as: 

G(E,M)p = — 7^G( E ,M)d 

2 1 

G(E,M)n = JjC(E,Af)d ~ ^G(E,M)b' (41) 

The up and down form factors, G^e,m)u and G(E,M)d 
are defined by convention to represent the up and down 
dressed quark form factors in the proton. The anomalous 
magnetic moments of the up and down quarks can be ex¬ 
pressed as k u = 2k p + k„, and Kd = n p + 2 n n , respectively. 
The Dirac and Pauli form factors for the up and down 
quarks can be written as: 

F(1,2)u = 2 F(i 2 )p + -E(l,2)rs 

F(l,2 )d = F(l,2 )p + 2-F( ll 2 )n- (42) 

The recent precision data on Gmu in the region of Q 2 
between 1.5 to 4.8 GeV 2 [78] and data on Geu/Gmu to 
Q 2 = 3.4 GeV 2 by Ref. [33] have enabled precise phe¬ 
nomenological fits to the proton and neutron form factors 
and detailed comparison to theory predictions. This allows 
one to extract information about the underlying contri¬ 
butions of the up and down quarks to the nucleon form 
factors. 

A calculation of the up and down quark form fac¬ 
tors from the available proton and neutron data was done 



Fig. 18. (Top panel) Up quark’s 1 Q 4 F 2u filled circle(black) 
and down quark’s K^ 1 Q 4 F 2 d filled triangle down (red) from 
Ref. 11911 . Up quark 1 Q 4 F 2u empty circle (black) and down 
quark’s K^ 4 Q 4 F 2 d empty triangle down (red) from Ref. 1192 :. 
(Bottom panel) Up quark’s Q 4 Fi u filled square (black) and 
down quark’s Q 4 F\d filled triangle (red) are from Ref. 1191 . Up 
quark’s Q 4 F\ U empty square (black) and down quark’s Q 4 Fid 
empty triangle (red) are from Ref. 11921 . . 


by Ref. m • The quark form factors were calculated 
to Q 2 = 3.4 GeV 2 by combining the measurements of 
Geu/Gmu by Ref. rrT31l3Tlf33inTTI[3?lH7(TI with the Kelly 
fit [53] to G Mn , Gm p and Ge p ■ Fig. [TH] is a plot of Q 4 Fi 
and Q 4 F 2 /nQ 2 versus Q 2 for the up and down quarks. 
The data is plotted at the Q 2 of the Geu/Gmu measure¬ 
ments and the error on the quark form factors is deter¬ 
mined by the error on the Geu/Gmu measurements. For 
Q 2 > 1.0 GeV 2 , the Q 2 dependence of both the F\ and F 2 
changes for the up and down quarks. For the up quark, 
Q 4 Ti and Q 4 F 2 /k continue to rise, while, for the down 
quark Q 4 Fi and Q 4 F 2 /k are plateauing or slightly drop¬ 
ping. 

Another separation of the up and down form factors 
was done by Ref. [192j . They used Ge p and Gm p from a 
extraction using cross section and polarization data which 
included two-photon exchange contributions [157] . In ad¬ 
dition, they added the data of Ref. [98] for Ge p /Gm p at 
low Q 2 . For the neutron form factors, they used the fit of 
Ref. [33] to Geu/Gmu and an updated parametrization 
of GMn using the data of Ref. [7^11 75 II76ll771178II 1 331 . The 
up and down form factors are calculated at the Q 2 of the 
proton data and are plotted in Fig. [TH] In the region of Q 2 
between 0.5 to 1.5 GeV 2 , F 2u from Ref. 11921 is slightly 
larger than that in Ref. [ 1911 and, correlated with that, 
F 2 d from Ref. )l92j is slightly smaller. In general, compar¬ 
isons between the two different flavor separations of the 
form factors give a sense that the size of the uncertainty 
due to two-photon exchange contributions and tensions in 
the data sets is relatively small and does not obscure the 
general trends in the Q 2 dependence of the up and down 
quarks form factors. 
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Great interest exists in extending the separation of 
quark flavors to higher Q 2 . In the spirit of Ref. m and 
Ref. [192j . we use the Kelly fit [55] for Gmu and Gm p 
while refitting HuGeu/Gmu and /. i p Ge p /Gm p to include 
the data since the Kelly fit was done. The fit form for 
l^nGEn/GMn IS- 

l^nGEn _ A\T _ /.q-v 

Gmu 1 + A 2 ^ + A 3 * T 

with A x = 2.6316, A 2 = 4.118 and A 3 = 0.29516. The fit 
form for fi p Ge p /Gm p is: 

H p Ge p _ 1 + B 0 t + Bit 2 + B 2 t 3 

Gm p 1 + B 3 t + B 4 T 2 + B 5 t 3 + Bqt a ’ 

with B 0 = -5.7891,= 14.493, B 2 = -3.5032, B 3 = 
-5.5839, B a = 12.909, B 5 = 0.88996 and B 6 = 1.5420. 
In Fig. [T9l the fits are compared to the world data for 
the proton and neutron F\ and F 2 /k. The shapes for the 
proton and neutron F 2 /k are nearly identical with the 
data on top on each other. A small shape difference in 
the proton and neutron F 2 /k dependence on Q 2 accounts 
for the difference between the up and down quark’s F 2 /k 
which is seen in Fig. [lH] 

Using the fits the nucleon form factors, the flavor sepa¬ 
ration can be extrapolated to higher Q 2 . To investigate the 



Q 2 (GeV 2 ) 

Fig. 19. Neutron and proton F\ and F 2 /k versus Q 2 . The data 
points are from the polarization experiments. The neutron F\ n 
and F 2n jKn data are plotted as filled diamond (black) and filled 
circle (red). The proton F\ v and F 2p /k p are plotted as filled 
square (green) and filled triangle (blue). The solid lines are 
calculations of Fi and F 2 based on the Kelly fit [55] for Gmu 
and Gm p and new fits to finGEn/GMn and h p Ge p /Gm p using 
Eqs. (ITU) and (ITU) . 


sensitivity of the extrapolation of the quark form factors, 
a different shape for Geu/Gmu can be used in calculating 
the quark form factors. As an example, the Geu/Gmu pre¬ 
diction from the Dyson Schwinger equation (DSE) model 
of Ref. 123] is plotted as a dash dotted line in Fig. [TT| 
The DSE models are discussed in the theory Sec. 14.21 
The up and down quark’s F\ and F 2 from the fit using 
Eqs. (PHI) and (1T51) are plotted in Fig. [2U] as a solid line in 
the region of Q 2 < 3.4 GeV 2 where Geu data exists, then 
extended as a dashed line when the fit is extrapolated 
to Q 2 = 12 GeV 2 . The fit shows that Q 4 F\d will have 
a zero crossing at Q 2 ~ 11.5 GeV 2 while Q 4 F Au starts 
to plateau above Q 2 of 7 GeV 2 . Both Q 4 F 2u and Q 4 F 2 d 
slowly drop-off above Q 2 of 3 GeV 2 with Q 4 F 2 d falling 
slightly faster. When the Geu/Gmu prediction from the 
DSE model is combined with the other form factors from 
the fit to calculate the quark form factor, then mainly the 
down quark’s form factors are modified. This is shown in 
Fig. [22] where the up and down quark’s F\ and F 2 using 
the Geu/Gmu prediction from the DSE model are plotted 
as a dash-dotted line. The zero crossing in Q 4 F\d moves to 
lower Q 2 . This demonstrates the need for precision mea¬ 
surements of all nucleon form factors to large Q 2 . Future 
experiments at JLab to extend the Q 2 range of the nucleon 
form factors measurements are discussed in Sec. [5] 



Fig. 20. Calculations of Fi and F 2 for the up and down quarks 
based on the Kelly fit [55] for Gm™ and Gm p and new fits to 
UnGEn/GMn and \i p Ge p IGm p using Eqs. flTU and QT]). The 
results from the fit are shown as a solid line in the region of 
Q 2 < 3.4 GeV 2 where Ge-h data exists, then extended as a 
dashed line when the fit is extrapolated to Q 2 = 12 GeV 2 . 
Calculations, which replace Geti/Gm-ti from fit using Eq. (1431) 
with Geu/Gmu from the DSE model of Ref. [123] while keeping 
the same Gmu, Gm p and [i p Ge p /Gm p fits, are plotted as a 
dash-dotted line. 


















V. Punjabi, C.F. Perdrisat, M.K. Jones, E.J. Brash, and C.E. Carlson: The Structure of the Nucleon 


23 


4 Theoretical Interpretations of Nucleon 
Form Factors 

We give an overview of theoretical work on nucleon elec¬ 
tromagnetic form factors. These form factors encode the 
information on the structure of a strongly interacting many- 
body system of quarks and gluons, such as the nucleon. 
This field has a long history and many theoretical at¬ 
tempts have been made to understand the nucleon form 
factors. This reflects the fact that a direct calculation of 
nucleon form factors from the underlying theory, Quantum 
Chromodynamics (QCD), is complicated as it requires, in 
the few GeV momentum transfer region, nonperturbative 
methods. Hence, in practice it involves approximations 
which often have a limited range of applicability. Despite 
their approximations and limitations, some of these non¬ 
perturbative methods do reveal interesting insight into the 
nucleon’s structure. 

Section 14.1.11 describes work simply fitting the form 
factors using an expansion in a transformed variable that 
allows convergence of a polynomial expansion for all values 
of Q 2 . These techniques have long been used in parame¬ 
terizing form factors measured in Weak decays, and part 
of the motivation here is to ensure a better extrapolation 
to obtain the slope at zero Q 2 , i.e., in determining the 
charge radius. 

Section 14.1.21 describes vector meson dominance of the 
photon-nucleon coupling, which in addition to being phys¬ 
ically motivated, provides good forms for fitting the nu¬ 
cleon form factors. Notable here is the fact that the falling 
G 1 ^ / G p m ratio was obtained from early on (although fits to 
the neutron form factors and form factors in the timelike 
region required tuning after data became available). 

Section 14.1.31 uses the ideas of dispersion relations to 
fit the form factors. In principle, if one knew the couplings 
and locations of all the poles and cuts in the NN chan¬ 
nels, one could calculate without further approximation 
the form factors at all Q 2 . In practice, the information is 
incomplete, and what one has are good parameterizations 
of the form factors that obey all the necessary analyticity 
properties and will converge everywhere. One then uses 
the available data to fit parameters in these functions, 
and hence obtain an accurate analytic representation of 
the data. 

Section f4.1.4l reviews the extensive work that has been 
done using constituent quark models to calculate the form 
factors. Relativity is crucial here, and many of the works 
use one of the Hamiltonian dynamical approaches enumer¬ 
ated by Dirac long ago 1193) . The constituent quarks are 
often thought to be representatives not only for elemen¬ 
tary quarks but also to represent not explicitly included 
contributions from gluons and higher Fock states, and as 
such may themselves have form factors that need to be 
parameterized. The outcome is usually a physically moti¬ 
vated form with parameters that need to be fit to data, 
and the fits in modern times are quite good for all the 
form factors. 

QCD has chiral symmetry, and chiral symmetry in our 
world leads to the existence of light mesons. These light 


mesons can then be part of the long range structure of 
the nucleons, and the pion cloud models that take these 
degrees freedom into account are described in Sec. 14.1.51 

To simply Fourier transform the electric and magnetic 
form factors to obtain the charge and magnetic densities 
is not valid relativistically, so this door to obtaining struc¬ 
ture information about the nucleons is closed. However, it 
can be shown that projections of the densities onto the 
transverse plane for a fast moving nucleon can be validly 
obtained from two-dimensional Fourier transforms of the 
Pauli and Dirac form factors. Work on transverse densities 
is described in Sec. 14.1.61 

A new approach to obtain approximate predictions for 
nucleon form factors is to use correspondences that have 
been discovered between gravitational theories in five di¬ 
mension and approximately conformal field theories, like 
QCD when considering only the light quarks. Results from 
this approach are described in Sec. 14.1.71 

The Dyson-Schwinger equations (DSE) are an infinite 
set of equations for the vertices, propagators, and other 
quantities related to calculating observables from a field 
theory. They can be truncated with some success, and 
they have been well studied in QCD, with DSE results for 
nucleon form factors presented in 14.21 

Perturbative QCD results for the form factors are on 
the face of it vitiated by the G^ / G P M measurements. The 
pQCD basics, and improvements to and extensions of the 
pQCD techniques, are reviewed in 14.3.11 

Form factors may also be obtained as integrals of gen¬ 
eralized parton distributions (GPDs), which are the am¬ 
plitudes for removing a quark from a nucleon and substi¬ 
tuting another quark with a different momentum and pos¬ 
sibly differing also in other quantum numbers. GPDs are 
important metrics of nucleon structure, and are measured 
in exclusive processes like 7 *N —> 'yN or q*7V —> pN. 
Their consequence for nucleon form factors is highlighted 
in Sec. 14.3.21 

Finally, nucleon form factor results from lattice gauge 
theory are reviewed in Sec. 14.41 The results so far are 
only for the isovector form factors, but the uncertainty 
limits are decreasing as further work is done, and lattice 
gauge theory has the advantages of being really QCD and 
not a model of QCD, and of obtaining results that are 
systematically improvable. 

4.1 Models of Nucleon Form Factors 

4.1.1 Conformal Fits to Form Factors 

There has been a lively discussion of what sorts of func¬ 
tions to use in fitting the form factors. Simple polynomial 
fits, for example, will not converge for moderate or large 
momentum transfers. The reason flows from the fact that 
the form factors are, from a mathematical viewpoint, an¬ 
alytic functions of their argument Q 2 = —q 2 1 except for 
cuts at known locations. The cuts are on the timelike side, 
and begin where one can have a photon to two pion tran¬ 
sition at q 2 = 4m^. A cut can be viewed as a weighted 
continuum of poles, so that there is a contribution to the 
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form factor containing a factor l/(q 2 — 4m 2 ). The weight¬ 
ing of this pole may be weak, but in principle its existence 
means that a polynomial expansion of the form factor will 
not converge for Q 2 > 4 m 2 . It is like the expansion of the 
geometric series 1/(1 — x), which does not converge for 

M > i- 

However, it is possible to make a mapping of Q 2 to an¬ 
other variable, denote it z, where a polynomial expansion 
in z is allowed. The trick is to find a transformation where 
spacelike momentum transfers all map onto the real line 
\z\ < 1 and timelike momentum transfers map onto the 
circle z = 1 (in the complex z-plane). Then since all poles 
of the form factors lie on the unit circle in z, a polynomial 
expansion in z is convergent everywhere inside the unit 
circle, i.e., for all spacelikc momentum transfers. 

This trick has been applied in the context of Weak 
interaction form factors, as for semileptonic meson decay, 
for some time (194111951 . It has now also been applied to 
fitting electromagnetic form factors by Hill and Paz 11621 
and by Lorenz et al. m- 

The variable z is given by the conformal mapping H62i 
1IM] . with t = q 2 , 


4.1.2 Vector Meson Dominance (VMD) 

The photon has the same J PC quantum numbers as the 
lowest lying vector mesons p(770), w(782), and </(1020). 
These mesons are prominent in the process e + e _ —> hadrons 
at the relevant timelike values of the CM energy squared 
q 2 > 0. One could hence expect that in elastic electron 
nucleon scattering at low spacelike momentum transfers 
q 2 < 0, some or much of the behavior of the coupling 
could be explained by the emitted photon converting to 
the strongly interacting meson which then attaches to the 
nucleon, as illustrated in Fig. [21] 


q 


/ P- w, cp 

ii 

- - - • - ► - 

Fig. 21. Vector meson dominance picture for the coupling of 
the photon (with four-momentum q) to a nucleon. 
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where Cut = 4m^ and one can easily enough verify that 
the mapping has the properties stated above. 

Fitting with a nonconvergent expansion can give good 
analytic fits to the data in any region where there is data 
to be fit. The danger lies in extending them outside the 
region where there is data. Such extrapolations can go 
awry, sometimes diverging wildly from physical expecta¬ 
tion and sometimes, depending on how far one extrapo¬ 
lates, there may be problems that are less visible. Here 
enters also the proton radius question, whose evaluation 
from form factors requires an extrapolation from the low¬ 
est Q 2 where there is data, down to Q 2 = 0. Extrapolating 
a fit made with an intrinsically convergent expansion is ar¬ 
guable safer. 

The two fits made to the electromagnetic form factors 
using the conformal variable z, however, differ in their 
conclusions regarding the proton charge radius. The ear¬ 
lier fit m used electron-proton scattering data avail¬ 
able before the Mainz experiment J59] published in 2010. 
They found a proton radius r p = 0.870 ± 0.023 ± 0.012, 
so their central value is closer to the CODATA value than 
to the muonic Lamb shift value. The other fit m used 
the 1422 data points from the Mainz experiment, and ob¬ 
tained Tp = 0.840 ± 0.015 fm (see also |196| ). 

Also in this section we may mention two recent reanal¬ 
yses of world e-p scattering data that obtain larger proton 
radii, one by Sick and Trautmann |1611ll97] and another 
by Graczyk and Juszczak m- The former was particu¬ 
larly mindful of effects of the charge density at large dis¬ 
tances upon the charge radius, and obtained r p = 0.886(8) 
fm, including all known data; the latter used a Bayesian 
framework and obtained r p = 0.879(7) fm, albeit without 
including the Mainz 2010 data m- 


That vector meson dominance (VMD) can explain much 
of the low Q 2 behavior of the form factors is confirmed by 
history. Before the p 109: . ui 2001 , and were ex¬ 

plicitly discovered in the early 1960’s, in reactions such as 
nN —> mrN or e + e _ —> pions, there were hints or pre¬ 
dictions of the existence gleaned from the behavior of the 
proton form factor. Nambu [202j in 1957 suggested that 
the observed form factor was consistent with the existence 
of a vector meson intermediary, and Frazer and Fulco [203] 
in 1959 in a famous dispersion analysis were more explicit, 
even suggesting a later confirmed mass range for the p me¬ 
son. 

A single vector meson exchange with simple couplings 
gives an my/{rriy — q 2 ) factor, from its propagator, for the 
falloff of the form factor. One can obtain a Q” 4 high mo¬ 
mentum falloff, in accord with observation or with pQCD, 
by having cancellations among two or more vector me¬ 
son exchanges with different masses, or more commonly 
in practice by giving the vector mesons themselves a form 
factor in their coupling to nucleons. 

An early example of a VMD fit to form factor data was 
given by Iachello, Jackson, and Lande |204l or IJL. They 
had several fits, but the one most cited is a 5 parameter fit 
with a more complicated p propagator that the form noted 
above, to account for the large decay width of the p meson. 
(The ui and </ are narrow enough that modifying their 
propagators gives no numerical advantage.) They in 1973 
predicted the falloff of Ge p /Ge p later seen experimentally, 
as illustrated along with some early data in Fig. [23] 

The IJL work was improved by Gari and Kriimpelmann 
[2()5H2061l207j to better match the power law pQCD expec¬ 
tations at high Q 2 , that F\ ~ Q -4 and F 2 ~ Q~ 6 , but also 
including some log Q 2 corrections to the falloffs based on 
the running behavior of the coupling a s (Q 2 ). 

Further improvement in VMD fits was made by Lomon 
[12611208] , who included a second p as the p'(1450), and 
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Fig. 22. Two early VMD fits are shown for nucleon form factor 
data from Iachello ef al. 12041 solid curve (blue), and Gari and 
Kriimpelmann (2051120611207! short dot-dashed curve (green). 
The data are from Hand et al. m open star (black), Price et 
al. [13] filled circle (magenta), Berger et al. [15] multiplication 
sign (red), and Bartel et al. [16] crossed diamond (cyan). 
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Fig. 23. Several VMD fits compared to the JLab Ge p /Gm p 
data. The solid curve (black) is the fit of Lomon 12081 . 
the dashed curve (blue) is that of Iachello, Jackson, and 
Lande 12041 . and the dotted curve (red) is that of Bijker and 
Iachello 12121 . 


later also a second u as the u/( 1419), and obtained a good 
parameterization for all the nucleon form factors. The first 
of the polarization transfer Ge p /Gmp measurements [7f 
was available in time for Lomon’s 2001 work i!26j . Lomon 
further tuned his fits pns| when the second set of polar¬ 
ization transfer data became available 0. 

Viewing the form factors as analytic functions of q 2 , 
the VMD forms are straightforward to analytically con¬ 
tinue to the timelike region (see, for example, Brodsky et 
al. nn or Dubnickova et al. [2T01 ). and compare to data 
that is now available. Workers in the field have done so, 
and have modified the VMD fits to give good accounts of 
data in both the timelike and spacelike region. 

In addition, the original IJL fits 12041 were not as good 
for the neutron as for the proton. Both the spacelike neu¬ 
tron form factors and timelike nucleon form factors were 
addressed in what may be termed IJL updates, by Iachello 
and Wan 21 lj and Bijker and Iachello [212j . both in 2004. 
Further, Lomon and Pacetti m have updated and ana¬ 
lytically continued the earlier Lomon fits in order to also 
give a good account of data in both timelike and spacelike 
regions. 

The VMD fits of course are fits to existing data, and 
they have been regularly updated as new data appeared. 
It will be interesting to check the “predictions” for the 
neutron form factors as newer data appears. A plot of the 
existing situation for protons is given in Fig. 1231 


4.1.3 Dispersion Analyses 

Newer works here include the dispersive analyses of the 
nucleon form factors by workers in Bonn [1640196111631 . 
There is also work that is only slightly older by Beluskin 
et al. and Baldini et al. [127LI214ll215j . The works include 
general analyses and fits to the form factors, as well as as¬ 
pects directly aimed at the resolution of the proton radius 
puzzle fTMimMUbSl . 

Dispersion relations relate the form factors in the space¬ 
like and timelike regions. More generally, the form factors 
are complex functions of q 2 that are analytic except for 
known cuts, and the form factors anywhere can be cal¬ 
culated if one knows just their imaginary parts at the 
cuts. The cuts are all on the real axis for timelike q 2 . 
The cuts run from q 2 ss 4to q 2 = oo. In practice, one 
cannot know the imaginary part of the form factors over 
this whole range, and uncertainty in knowing the form 
factors in the timelike region builds nonlinearly to larger 
uncertainty in predicting the form factors in the spacelike 
region, especially at higher Q 2 . 

At lower Q 2 , one specific boon of the dispersive treat¬ 
ment is that the connection between the timelike and 
spacelike regions puts an extra constraint on the form fac¬ 
tors and their slope at spacelike threshold. This means 
that the determination of the charge and magnetic radii 
are not purely extrapolations of the scattering data, but 
is an interpolation in this procedure and hence arguably 
more reliable. 
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A technique involving dispersion relations, used by both 
m and by is to parameterize the imaginary part 

of the timelike form factors, and determine the parame¬ 
ters by making a least squares fit to the known spacelike 
and timelike data. One thereby obtains a representation 
of the form factors that one can use in regions where there 
is not yet data. 

Ref. [2141 . published in 2006, applies dispersion re¬ 
lations to the ratio Ge p /Gmp , with the assumption of 
no zeros in Gm p ■ They used some large-uncertainty-limit 
Gep/Gmp data in the timelike region, obtained from an¬ 
gular distributions in e + e _ —> pp or the reverse, to sup¬ 
plement the polarization data in the spacelike region. One 
of the main goals was to compare to models that fit the 
spacelike data, especially to the continuations of those 
models to the timelike regions HE]. 

They found that there was a zero in Ge p {i ? 2 ) at about 
11 GeV 2 spacelike momentum transfer squared, and found 
that the Phragmen-Lindelof theorem, which leads to the 
statement that \GE P {q 2 ) / GM P (q 2 )\ should be the same at 
very large momentum transfers, whether spacelike or time¬ 
like, was satisfied, albeit with opposite signs. Many of the 
more purely phenomenological models differed on the lat¬ 
ter point. This work m preceded the completion of the 
polarization experiments at the highest current Q 2 ms 
nun . The dispersive aspects have not been updated in 
more recent works by some of the same authors, e.g. EE], 
but one can see that the results would not be materially 
changed by the newer data. 

Lorenz et al. ira in their update of 027, apply the 
dispersion analysis to Ge p and Gm p separately. An im¬ 
portant improvement in the newer work m is that it 
includes the recent Mainz data [5Uj in its fit. A salient 
outcome of this analysis is that the proton charge radius 
comes out at a value r p = 0.84 (1) fm, in agreement with 
the value found in the muon hydrogen Lamb shift mea¬ 
surement. 


4.1.4 Constituent Quark Models 

Constituent quark models (CQMs) have been used to un¬ 
derstand the structure of baryons, beginning when quarks 
were first conjectured and predating the establishment of 
QCD as the theory of the strong interactions. Indeed, the 
observed spectroscopy, particularly the existence of the 
Z\ ++ , played a crucial role in bringing to light the quan¬ 
tum number of color. In the CQM, the nucleon is a quan¬ 
tum mechanical ground state of three quarks in a con¬ 
fining potential. More generally, ground state baryons are 
composed of three quarks, selected from up (u), down (d) 
and strange (s) flavors, and are described using 517(6) 
spin-flavor wave functions and a completely antisymmet¬ 
ric color wave function. 

Early CQMs concentrated on explaining static prop¬ 
erties, including magnetic moments and transition ampli¬ 
tudes. Examples are the models of De Rujula, Georgi, and 
Glashow m and of Isgur and Karl m■ In the latter, 
the quarks were in a harmonic oscillator potential, and 
at least at first the wave functions were nonrelativistic 


product wave functions, and the ground state baryons ap¬ 
peared as a 56-plet of 517(6). QCD by this point having 
been discovered, the hyperfine splittings, e.g., between the 
nucleon and the Z\(1232), were given by a one-gluon ex¬ 
change potential added to the confining potential. The 
one-gluon exchange also generates a small tensor interac¬ 
tion that leads to some D-state admixture into the ground 
state baryons. This in turn allows some non-zero electric 
quadrupole (E 2) and Coulomb quadrupole (C2) nucleon 
to Z\(1232) transitions, in accord with observation. 

However, form factors require a relativistic treatment. 
At high Q 2 , nonrelativistic treatments lead to form factors 
that are far too small compared to data. At low Q 2 , the 
charge radius defined from the slope of the form factor 
has contributions corresponding to the RMS charge radius 
known from nonrelativistic treatments, but also has recoil 
terms proportional to the Compton wavelength or inverse 
mass, squared, of the target. The latter are absent in any 
nonrelativistic model, and the nucleon is light enough for 
this to be a problem. 

A crucial question for a form factor calculation, since 
the nucleon must be moving after or before the interac¬ 
tion or both, is how the wave function in the rest frame 
transforms to a moving frame. This is not a trivial ques¬ 
tion, and the answer can, depending on the formalism, 
be dependent on the interactions binding the quarks. For¬ 
mally, one needs to know how the eigenfunctions of the 
mass and spin operators can be viewed as unitary repre¬ 
sentations of the Poincare group, whence it will be known 
how they change under Poincare, which includes Lorentz, 
transformations. The generic ways this can be done were 
laid out by Dirac 1193] . There are three forms of dynam¬ 
ics, which are the instant, point, and light-front forms. 
These differ in which generators form the kinematic sub¬ 
group of the Poincare group. This is the subgroup whose 
commutators do not involve the interactions among the 
constituents. The Poincare group has ten generators, four 
space-time translations (momentum operators), three spa¬ 
tial rotations, and three boosts. In a given representation 
each of these may be kinematical or interaction depen¬ 
dent, or dynamical, meaning dependent on the specifics 
of the interaction potential. The latter cannot usually be 
dealt with in an exact way, but must be dealt with ap¬ 
proximately or numerically in a practical calculation. 

The point form has all boosts and rotations kinemat¬ 
ical, meaning that as operators in a field theory they can 
be written out without having to know the interaction 
Lagrangian or interaction Hamiltonian. Straightforwardly, 
the angular momenta and Lorentz boost are the same as 
in the free case. However, all four components of the mo¬ 
mentum operators are interaction dependent in this case. 

The instant form has the rotation operators and space 
components of the momenta kinematical. Eigenstates of 
the angular momentum are then easy to construct. How¬ 
ever, the time component of the momentum, or Hamilto¬ 
nian, and the boosts are dynamical. Boosts, then, require 
knowing and including effects of the interaction in order 
to ascertain important infomation, e.g., the momentum 
space wave function, of the boosted state. 
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The light-front form has seven kinematical generators. 
This is the maximum possible. The three dynamical gener¬ 
ators are one component of the four-momentum operator 
(for which the mass operator obtained from M 2 = p 2 ^ may 
be substituted) and two light-front transverse rotations 
(or light-front transverse boosts, meaning here two par¬ 
ticular linear combinations of the two ordinary transverse 
rotations and the two ordinary transverse boosts). Light- 
front calculations have the advantage that states can be 
easily and exactly transformed from one frame to another 
the using the (kinematic) longitudinal boost, and the two 
kinematic light-front transverse rotations (the linear com¬ 
binations orthogonal to the one previously mentioned). 
However, light-front calculations have the disadvantage 
that states of definite angular momentum are difficult to 
construct because the rotation operators are interaction 
dependent. 

A calculation of the form factors also requires know¬ 
ing the electromagnetic current operator. It is usually as¬ 
sumed that the photon only interacts with one quark in 
the nucleon. 

The relative ease of exactly transforming states from 
the frame where the wave functions are calculated or oth¬ 
erwise given, to any other frame, makes the light-front 
form attractive for form factor calculations. The light- 
front form in this context was introduced by Berestetsky 
and Terentev [218H219J , and later developed by Chung and 
Co ester I 22ffl . In these calculations one begins with some 
wave function that has been developed in CQMs designed 
to study the baryon spectrum. The light-front form of the 
wave function is obtained by a Melosh or Wigner rotation 
of the Dirac spinors for each quark. These come about 
because the usual CQM models use spinors that in mo¬ 
mentum space are obtained from rest spinors by a direct 
boost, while the light-front spinors are obtained by a longi¬ 
tudinal boost followed by a kinematical light-front trans¬ 
verse boost. Undoing one and then applying the others 
leads to the same momentum, but leaves a residual rota¬ 
tion that does not affect the momentum, but does affect 
other quantities such as the spin. While the Melosh rota¬ 
tion is not difficult conceptually, the expressions it leads 
to are tedious to write out. 

If, in addition, one calculates in the so-called light- 
front (or Drell-Yan) frame characterized by q + = 0, then 
momentum conservation ensures that the current matrix 
elements connect only states whose Fock components have 
the same number of constituents. There are, for example, 
no matrix elements connecting qqq to q 4 q states, so that a 
consistent calculation can be done using only three-quark 
states. 

Chung and Co ester used Gaussian wave func¬ 

tions. They did obtain a falling Ge p /Gmp ratio. This ap¬ 
parently m is a feature shared by many relativistic cal¬ 
culations and occasioned by the Melosh transformation. 
However, the form factors fell far too fast at large Q 2 . 
Schlumpf [22111222] used instead a wave function with a 
power law falloff, fitting parameters in his wave function 
to static baryon properties. The high Q 2 falloff was now 
in line with data, including some at that time new neu¬ 


tron data 122: and also showed a i^p/Tip ratio that fell 
more slowly than 1 /Q 2 1 in qualitative agreement with later 
data. 

French, Jennings, and Miller [22311224] focused on the 
effects of the nuclear medium upon the form factors, but 
also calculated the form factors of single free nucleons. 
They followed the work of Schumpf [22111222] and simi¬ 
larly found a decreasing Ge p /Gm p ratio, obtaining a zero 
between 5 and 6 GeV 2 for Q 2 . 

Improvements in the detailed quality of the fit can 
come by introducing Dirac and Pauli form factors for the 
quarks. This can be justified by arguing that nucleons are 
not just bound states of three quarks, but have further 
constituents in the form of gluons and quark-antiquark 
pairs. Modeling the nucleon with three quarks means the 
constituent quarks are also representing the additional 
components present in a complete Fock space expansion 
of the nucleon, and this gives them an effective structure 
represented by the quark form factors. This viewpoint was 
taken by Carderelli et al. [22511226] (plus further references 
contained in the latter) to produce good fits to both the 
proton and neutron form factor data then available. They 
used the light-front formalism and quark wave functions 
obtained from a potential of Capstick and Isgur [ 2271 , and 
made the point that the one-gluon exchange is crucial to 
obtaining sufficient high momentum components in the 
wave function to explain the form factor data. 

A different starting wave function, now in the context 
of the point form formalism, appears in the hypercentral 
constituent quark model of De Sanctis, Santopinto, and 
others ( [22811225] and references therein). The feature here 
is that the confining potential is treated as a function of an 
average separation defined by the RMS sum of the quark 
positions, relative to the CM, and there is also a term to 
give the hyperfine splitting. Parameters of the potential 
are fit to the baryon mass spectrum. With the inclusion 
here also of form factors for the constituent quarks, good 
fits are obtained for the nucleon form factors, with up¬ 
dates 12291 succeeding the latest polarization transfer Gl 
results [IB]. 

A comparable amount of high-momentum components 
in the nucleon wave function was obtained in the Goldstone- 
boson-exchange (GBE) quark model [23011231] . This model 
relies on constituent quarks and Goldstone bosons, which 
arise as effective degrees of freedom of low-energy QCD 
from the spontaneous breaking of the chiral symmetry. 
The resulting CQM assumes a linear confinement poten¬ 
tial supplemented by a quarkquark interaction based on 
the exchange of pseudoscalar Goldstone bosons, which is 
the source of the hyperfine interaction. It was shown I23U1 
I25T1 that the GBE CQM yields a unified description of 
light- and strange-baryon spectra. The GBE CQM was 
used in [23211233] to calculate the nucleon e.m. form factors 
in the point-form. The neutron charge radius is well de¬ 
scribed in this model and is driven by the mixed-symmetry 
component in the neutron wave function. In contrast to 
the light-front calculation 122611234] . it was found that 
when performing a point-form calculation of the nucleon 
e.m. form factors at larger Q 2 within the impulse ap- 
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proximation, i.e. considering only single-quark currents, 
a surprisingly good overall description of the nucleon e.m. 
form factors can be obtained, using point-like constituent 
quarks only. When looking at details of Refs. [252112331 . 
the agreement is worse though for Gm p which is under¬ 
predicted at larger Q 2 , and the ratio of Ge p /Gmp is over¬ 
predicted at larger Q 2 , see Fig. l24l Similar findings have 
also been obtained in the point-form calculation of [235 1 
for the one gluon exchange CQM. The overall success of 
the point-form result using point-like constituent quarks 
was attributed in |232ll233ll235] to the major role played 
by relativity. Such a finding is remarkable in view of the 
expected finite size of the constituent quarks, as discussed 
above. 

An explanation for the above finding for the nucleon 
e.m. form factors in the point form, using the single-quark 
current approximation, has been suggested by Coester and 
Riska I23H| . When the spatial extent of the three-quark 
wave function is scaled (unitarily) to zero, both instant 
and front forms yield form factors independent of the mo¬ 
mentum transfer. Therefore, to reproduce the experimen¬ 
tal fall-off of the nucleon e.m. form factors at large mo¬ 
mentum transfers requires the introduction of constituent 
quark form factors. In contrast, when the wave function 
in point form is scaled unitarily to zero (so-called point 
limit), a non-trivial scaling limit is obtained for the form 
factors, depending on the shape of the wave function. At 
high values of momentum transfer, the scaled form factors 
decrease with an inverse power of the momentum transfer. 
The power is determined by the current operator and is 
independent of the wave function. An explicit compara¬ 
tive calculation of the baryon e.m. form factors between 
the three different forms was performed in 12371 using a 
simple algebraic form for the three-quark wave function, 
depending on two parameters. It was verified that a qual¬ 
itative description of the nucleon form factors data de¬ 
mands a spatially extended wave function in the instant- 
and front-form descriptions, in contrast to the point-form 
description which demands a much more compact wave 
function. 

A manifestly covariant CQM calculation within the 
Bethe Salpeter formalism and using an instanton-induced 
interaction between quarks has been performed by Merten 
et al. 238] . Although this model reproduces the baryon 
spectrum, it can only qualitatively account for the Q 2 de¬ 
pendence of the nucleon e.m. form factors. 

Another type of covariant CQM calculation was done 
by Gross, Ramalho, and Pena 12331 . partly based on ear¬ 
lier work of Gross and Agbakpe avoiding questions 
of dynamical forms by staying in momentum space. They 
use a covariant spectator model, where the photon inter¬ 
acts with one quark and the other two quarks are treated 
as an on-shell diquark with a definite mass. They too take 
the view that the constituent quark includes, at least at 
lower Q 2 , effects from higher Fock states and so the quark 
itself should have a form factor, including the possibility of 
nonzero quark anomalous magnetic moments. They note, 
as others have, that obtaining a good fit to the neutron 
electric form factor Geu requires isospin breaking. Since 


they use VMD forms for the quark form factors, the dif¬ 
ference between the p (isospin-1) and the uj (isospin-0) 
couplings is one sufficient way to obtain this. They obtain 
forms with parameters that they fit to the data. Their fit 
from the 9-parameter “model IV” is quite precise. 

Of interest, especially noted in | i239l although perhaps 
somewhat relaxed in m , also by Gross et al ., is that only 
5-state quark wave functions are needed to obtain good 
fits to the form factors. Other work (the article of Brod¬ 
sky and Drell (242] is an early example) argues, on the 
other hand, that nonzero angular momenta are necessary 
for fitting the form factor data. This can be a difference 
in organization of the calculation. Looking in particular 
at the magnetic moments, if there are only 5-states, the 
anomalous magnetic moments of the nucleons can only 
come from intrinsic anomalous magnetic moments of the 
quarks, which can be nonzero if quark form factors are 
allowed; whereas if the quarks are all elementary par¬ 
ticles their anomalous magnetic moments are zero, and 
higher Fock states or nonzero orbital angular momenta 
are needed to produce the observed magnetic moments of 
the nucleons. 
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Fig. 24. The JLab Ge p /Gm p data compared to the results 
of a selection of constituent quark models. The short dashed 
curve (blue) is from Boffi et al. 12331 . the solid (orange) from 
de Melo et al. 12431 . the long dash (magenta) from Gross et 
al. 1239] . the dotted (red) from Chung and Coester [220], and 
the dash-dot (cyan) from Cardarelli et al. [225j . 



4.1.5 Pion Cloud Models 

In nature, the up and down quarks are nearly massless. 
In the exact massless limit the QCD Lagrangian is sepa- 
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rately invariant under rotations in flavor space of the left 
handed and right handed quarks. I.e., there is invariance 
under a chiral symmetry group SU(2)l x SU(2)r. Parity 
is also a symmetry of the Lagrangian, and fully compat¬ 
ible with the chiral symmetry, and so it would seem that 
for any state of QCD, there would be an equal mass state 
of opposite parity. However, parity doubling is not seen 
in nature. Instead we find the alternative, which is that 
there exist massless negative parity scalars, the Goldstone 
bosons. For two quark flavors, there are three Goldstone 
bosons, the pions, which acquire a mass in nature because 
of the symmetry breaking due to the quark masses. Since 
the pions are light, they dominate the long distance behav¬ 
ior of the nucleon wave functions, and have a potentially 
important effect on the low-momentum transfer behav¬ 
ior of hadronic form factors. Hence one can improve the 
constituent quark models by including pionic degrees of 
freedom. 

One early improvement to constituent quark models 
was in the context of the bag model of hadrons, where 
a number of workers, in or about 1980, including Brown 
and Rho [244] , Jaffe [245j . and Thomas, Tlieberge, and 
Miller [2461 secured the chiral symmetry of the model by 
including coupling to pions. The vision was then of quarks 
within the boundary of the bag and pions without, and the 
model was called the Cloudy Bag Model. However, states 
of the bag model are expressed in coordinate space as set 
of independent particle wave functions for each quark. 
Turning these states into momentum eigenstates is not 
a trivial problem in general, because there are center-of- 
rnass fluctuations that must be removed, and momentum 
eigenstates are needed to discuss the form factors. Addi¬ 
tionally, the bag has a simple spherical boundary in the 
rest frame, and it must be understood how the states alter 
under Lorentz transformation in order to make accurate 
form factor calculations for nucleons. Hence some time 
passed before Lu et al. in 1998 [24T calculated nucleon 
electromagnetic form factors in the cloudy bag model. 
They used the Pierels-Thouless formalism and a plausi¬ 
ble hypothesis for the effects of Lorentz transformations 
(details in [2171 ). and obtained a good description of the 
nucleon electromagnetic form factors for Q 2 < 1 GeV 2 . 

Miller extended the calculation to larger Q 2 using a 
light-front version of the cloudy bag model calculation 
m- Starting from constituent quarks [224], using the 
Schlumpf wave function instead of bag wave functions for 
the quark core, Miller calculated the effects of the pion 
cloud through one-loop diagrams. The model gives a rela¬ 
tively good overall account of the form factor data at both 
lower Q 2 and higher Q 2 . 

The cloudy bag model is one chiral quark model which 
treats the effect of pions perturbatively. Other quark mod¬ 
els which calculated nucleon electromagnetic form factors 
using perturbative pions can be found e.g. in the early 
works of |248U249H250j , as well as in the already discussed 
works of [23011231] . Recently, the above chiral quark mod¬ 
els where pions are included perturbatively have been im¬ 
proved in [2511 . This work extends a previous work of 
Esa by dynamically dressing bare constituent quarks by 


mesons to fourth order within a manifestly Lorentz co¬ 
variant formalism. Once the nucleon and A hyperon mag¬ 
netic moments are fitted, other e.m. properties, such as 
the nucleon e.m. form factors at low momentum transfers, 
follow as a prediction. It was found in [2511 that the me¬ 
son cloud is able to nicely describe the form factor data in 
the momentum transfer region up to about 0.5 GeV 2 . To 
extend the calculations to larger Q 2 , a phenomenological 
approach has been adopted in [251] by introducing bare 
constituent quark form factors which were parameterized 
in terms of 10 parameters. Such parameterization makes it 
plausible to simultaneously explain the underlying dipole 
structure in the nucleon e.m. form factors as well as the 
meson cloud contribution at low Q 2 which results from the 
underlying chiral dynamics. In a later paper , a model 
calculation for the bare constituent quark form factors has 
been performed and applied to the electromagnetic prop¬ 
erties of the TV —» A transition, and also of the TV —> Roper 
transition [254] . The latter paper also includes updates of 
that group’s fits to the nucleon form factors. 

When pion effects dominate nucleon structure, their 
effects have to be treated non-perturbatively. A nonper- 
turbative approach which has both quark and pion de¬ 
grees of freedom and interpolates between a CQM and the 
Skyrrne model (where the nucleon appears as a soliton so¬ 
lution of an effective nonlinear pion field theory) is the 
chiral quark soliton model (xQSM). As for the Skyrrne 
model, the ,\QSM is based on a 1/TV C expansion (with 
N c the number of colors in QCD). Its effective chiral ac¬ 
tion has been derived from the instanton model of the 
QCD vacuum [2551 . which provides a natural mechanism 
of chiral symmetry breaking and enables one to gener¬ 
ate dynamically the constituent quark mass. Although in 
reality the number of colors TV C is equal to three, the ex¬ 
treme limit of large TV C is known to yield useful insights. 
At large N c the nucleon is heavy and can be viewed as TV C 
“valence” quarks bound by a self-consistent pion field (the 
“soliton”) 12561 . A successful description of static proper¬ 
ties of baryons, such as mass splittings, axial constants, 
magnetic moments, form factors, has been achieved (typ¬ 
ically at the 30 % level or better, see [ 2571 for a review 
of early results). After reproducing masses and decay con¬ 
stants in the mesonic sector, the only free parameter left 
to be fixed in the baryonic sector is the constituent quark 
mass. When taking rotational (1/TV C ) corrections into ac¬ 
count, this model achieved a qualitative good description 
of the nucleon electromagnetic form factors in the range 
Q 2 < 1 GeV 2 , using a constituent quark mass around 
420 MeV [258] . The chiral soliton model obtained a de¬ 
crease of the Ge p /Gm p ratio with increasing Q 2 already 
in the late 1990’s. 

Holzwarth [259] extended the chiral soliton model by 
including the p and to meson propagators for the isovec¬ 
tor and isoscalar channels, respectively. Furthermore, to 
extend the range in Q 2 of the predictions, he uses a rela¬ 
tivistic prescription to boost the soliton rest frame densi¬ 
ties to the Breit frame. Such prescription is also used to 
extract radial charge and magnetization rest frame densi¬ 
ties from experimental form factors, as will be discussed 

































30 


V. Punjabi, C.F. Perdrisat, M.K. Jones, E.J. Brash, and C.E. Carlson: The Structure of the Nucleon 


in Sect. 14.1.61 Using 4 fit parameters (one effective boost 
mass and three free parameters to fix the couplings of p 
and u> mesons), the model was found to provide a good 
account of the detailed structure of the nucleon e.m. form 
factors in the low Q 2 region. In particular, for Ge p /Gmp 
it predicts a decreasing ratio in good agreement with the 
data. At larger Q 2 , the boost prescription gives a rea¬ 
sonably good account of the data (except for Gmu ) and 
predicts a zero in Ge p around 10 GeV 2 . Due to the uncer¬ 
tainty introduced from the particular choice for the boost 
prescription, the high Q 2 behavior (for Q 2 larger than 
about 4m 2 ) of the e.m. form factors is however not a pro¬ 
found prediction of the low-energy effective model. 

Cloet and Miller [260] . in addition to fitting the elec¬ 
tromagnetic form factors with a quark plus pion cloud 
model, have the further goal to accurately reproduce the 
spin fraction of the proton that comes from quark spin. 
It will be remembered that the EMC collaboration found 
that little of the proton spin came from quark spin, and a 
more modern analysis [261] gives the quark spin fraction 
of the proton spin as 36.6jj^g% (for cc m i n = 0.001 and 
la). The Cloet-Miller model uses a light-front formalism 
with a quark-diquark system accounting for 71% of the 
nucleon state, by probability, plus a quark-diquark core 
with a pion accounting for the rest. Gluons are not ex¬ 
plicitly included. They have 10 parameters, which they fit 
to the electromagnetic form factor data, obtaining a good 
representation of the data and/or the Kelly fit thereto. 
They have a zero of Ge p at Q 2 = 12.3 GeV 2 . The to¬ 
tal quark plus diquark spin they obtain is 36.5% of the 
proton’s h/ 2, in fine accord with expectation. 

4.1.6 Transverse Densities 

Nonrelativistically, form factors and charge distributions 
are Fourier transforms of each other. Relativistically, this 
is no longer true, because of quandaries in transforming 
the nucleon wave function from one reference frame to an¬ 
other. In terms of the charge or magnetic radius this leads 
to additive terms that may be called recoil terms, pro¬ 
portional to Q 2 /m 2 N with coefficients that are not trivial 
to calculate. For heavy targets, these terms may be ig¬ 
nored, but the nucleon is light enough that precision work 
cannot proceed nonrelativistically. One cannot obtain ac¬ 
curate spatial images of the nucleon charge or magnetic 
densities by just Fourier transforming the charge or mag¬ 
netic form factors. 

However, one can obtain kinematically correct, and ac¬ 
curate to the same level that the input data is accurate, 
charge distributions if one is willing to adopt a new view¬ 
point. The viewpoint is that of a light front moving to¬ 
wards a nucleon, or equivalently of an observer viewing 
a nucleon approaching at nearly the speed of light. The 
charge density that will be seen is two dimensional, or the 
3D charge density of the nucleon projected onto a plane 
transverse to its direction of approach. Additionally, the 
charge density that will be seen is not the charge density 
obtained from the electromagnetic current component J°, 
but rather densities coming from J + = J° + J 3 . 


Typically one chooses the z-direction as special, choos¬ 
ing it along the direction P = {p+p')/ 2, where p and p' ate 
the incoming and outgoing nucleon momenta, and further 
arranges the frame so the photon momentum has q + = 0, 
and its transverse (lying in the xy plane) momentum is 
denoted q±. The charge density projected onto the trans¬ 
verse plane is shown in works of Burkardt and of Miller 
and of others |2621l2631l2641l265j to be 

x(P+, ^,A| J+( 0 )|P+,-^,A), (46) 

where b is the position in the xy plane relative to the 
nucleon’s CM, A is the (light-front) helicity, and the lon¬ 
gitudinal and transverse components of the incoming and 
outgoing nucleon’s momenta are indicated separately. For 
the density pqn of an unpolarized nucleon one finds, fol¬ 
lowing Miller [2631, that 

P0N(h) = J™ ^QM b Q) F i(Q 2 ), (47) 

where Q 2 = qy, Jo is the Bessel function, and Jj is the 
Dirac form factor. Further, if one polarizes the nucleon 
transversely in the ^-direction, one obtains EMI, 

PTN(b) = P0N{b) 

<48> 

where 4>b is the azimuthal angle of b , and F -2 is the Pauli 
form factor. 

The result for the unpolarized density distribution of 
the proton is not startling and is not shown here, but the 
result for the neutron is quite striking 13631 . The neutron 
charge density is found to be negative near its center as 
shown in Fig. 1251 it had long been known to be negative 
far from the center. 

2nbp Qn [1/fm] 



Fig. 25. The 2D projected charge density of a neutron mov¬ 
ing rapidly towards an observer, shown as a function of the 
distance from the CM of the neutron. 
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For a transversely polarized neutron, one finds a charge 
separation as shown in Fig. | 


based on figures in 
The upper panel shows the difference between the charge 
densities of the polarized and unpolarized neutron. One 
sees negative charge above and positive charge below. The 
lower panel gives similar information by comparing the po¬ 
larized and unpolarized charge distributions along a single 
line, the y-axis. (One should know that time reversal in¬ 
variance forbids an electric dipole moment for a stationary 
elementary particle, but for a moving particle the electric 
dipole moment d is given by d = v x fi, where v is the 
velocity of the particle and fi is its magnetic moment.) 


field theory (AdS/CFT) correspondence, initiated by Mai 
dacena 
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Fig. 26. Upper panel: the 2D projected charge density pj f° r 
a neutron polarized in the ^-direction (to the right, in the fig¬ 
ure). Dark areas represent negative charge, light areas, positive 
charge. Lower panel: the charge densities along the y-axis for 
a neutron polarized along the ^-direction; po„ is given by the 
dashed curve (blue), and pm is given by the solid curve (red). 


where the most relevant approximately con¬ 
formal field theory is QCD and we talk of the AdS/QCD 
correspondence, where early applications were made by 
Erlich et al. [2B71 . The idea is that some string theory in 
ten dimensional space has a symmetry that is divided so 
that for five dimensions we are and remain in the ground 
state in those dimensions, and the theory in the other 
five dimensions mimics a gravitational theory in a five 
dimensional anti-de Sitter space, or AdSs. The AdSs pos¬ 
sesses a SO( 2,4) symmetry, which is important because it 
is the same symmetry possessed by a conformal (in prac¬ 
tice, a theory with all masses zero) field theory in four di¬ 
mensions. One can exploit the sameness of the symmetry 
group to find numerically viable relations between the so¬ 
lutions to the gravitational theory in five-dimensions and 
the conformal field theory in four-dimensions. For an ex¬ 
tensive review, see [2681 . 

Since QCD is not a conformal theory, maintaining a 
correspondence requires also breaking the symmetry of 
the AdS space, which in practice is done either with the 
hard wall model, where the AdS space is cut off at long 
distances in the fifth dimension, or by the soft wall model, 
where an extra potential is introduced that suppresses 
long distance propagation in the fifth dimension. 

Much of the work on individual particles has focused 
on the bosonic sector, studying both quarkic hadrons and 
glueballs of various spins, and obtaining masses, decay 
constants, and charge radii. The actual correspondence 
is between operators in the four-dimensional space, such 
as the electromagnetic current or the energy momentum 
tensor T^, and fields in the five-dimensional space with 
corresponding quantum numbers, such as a vector field or 
the metric g iJV . Results can be compared to experiment at 
the few times 10% level (268 ;. 

Obtaining results for fermions is more involved than 
for bosons. One approach is to build from the bosonic 
sector, and treat the fermions as Skyrmions within that 
model [26011270] . Another approach is to begin in AdS with 
fundamental fermions that interact with the AdS gravita¬ 
tional background [27I11272y2781l274y276112761l2771[mi . In 
this version the five-dimensional Lagrangian is still rela¬ 
tively simple. It has terms for gravity (the scalar curvature 
and the cosmological constant terms) and for vector fields, 
and the terms for the fermions are simple given the con¬ 
text of interactions with gravity and the vector field, 


C F = y/g(^e%r A D N &-^(D N &yr°e%r A &-MMy 

( 49 ) 

for the hard wall version of the model, and the covariant 
derivative is Djv = 8n + ^nab[B a , r B ] — iVN , where Vn 
is the vector field which will be dual to the electromagnetic 
current when we make the correspondence. The indices N, 
A, B , etc., each run over five values as appropriate to the 
five-dimensional space. In the above Lagrangian, gravity 
enters via the metric <7 mjv and is seen in its determinant 
g and also in the spin connection lonab , which will not be 
mate QCD results is to use the anti-de Sitter space/conformalgiven in detail here. 


4.1.7 Correspondences with Higher Dimensional Theories 

A recent, exotic, and interesting way to obtain approxi- 
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As an example, results obtained using the AdS/CFT space as 
correspondence for the proton electromagnetic form fac¬ 
tors by Abidin et al. 1276j are shown in Fig. [27] for both 
the hard wall and soft wall models. 


S F (p 2 ) = 


iF(p 2 ) 

$ + M(p 2 ) ’ 


(50) 




Fig. 27. Results for the proton electromagnetic form factors 
obtained using the AdS/CFT correspondence [276] . with the 
soft wall model in dashes (orange) and the hard wall version in 
dot-dash (violet), compared to the fits of Arrington et al. 12791 
shown as solid lines (blue). 


4.2 Dyson-Schwinger Equations and Diquark Models 

The Dyson Schwinger equations (DSE) are generically a 
non-perturbative approximation for obtaining results for 
a field theory, in the present case QCD. The equations are, 
in principle, an infinite set of coupled integral equations. 
In practice, they must be truncated, in a way that pre¬ 
serves all symmetries of QCD, in order to proceed with 
any calculation. For a general DSE review, see 12801 . 

One accomplishment of the DSE follows the solution 
for the full quark propagator, represented in momentum 


where the normalization F{p 2 ) and the mass M{p 2 ) be¬ 
come functions of momentum because of interactions. With 
relatively simple truncations and modeling of the QCD in¬ 
teractions, the DSE obtain a mass function in good agree¬ 
ment with lattice calculations. The mass M(p 2 ) is several 
hundred MeV at small p 2 and falls smoothly to the small 
values at large p 2 that one might expect in perturbation 
theory. 

Also early in the DSE program is building a model 
of the quark-quark and quark-antiquark interactions that 
will reproduce data on, among other quantities, the pion 
mass and decay constant. 

One then uses the same quark-quark interactions de¬ 
veloped in meson studies to obtain a three quark wave 
function model for the nucleon by solving the three-cody 
Fadeev equations. There arise significant diquark contri¬ 
butions, i.e., significant quark-quark correlations, which 
have a strong effect on the form factors one obtains. Since 
the quarks in this model are dressed, many of their fea¬ 
tures are different from expectations for pointlike fermions. 
One finds in particular large quark anomalous cliromo- 
magnetic moments, which affect the quark-gluon interac¬ 
tions, which lead to large quark anomalous magnetics mo¬ 
ments in the quark-photon interactions, which in turn are 
needed to obtain good fits to the nucleon electromagnetic 
form factor data. 

The theoretical DSE results, from the work of Cloet, 
Roberts, and others )281ll282[l283j . show a falloff of the 
G p e /G p m ratio similar to what is seen in the data; see 

Fig.m 
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Fig. 28. An illustration of what may be learned about bound 
state substructure, in the context of a DSE model 12821 . from 
the measured Ge p /Gm p ratio. Parameter a measures how 
quickly the dressed or constituent quark mass approaches its 
asymptotic or perturbative value; a. = 1 is a benchmark favored 
by fits to meson masses, and larger values of a accelerate the 
rate of approach to the asymptotic quark masses. The data is 
from [gEmm!. 
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Qualitatively, the behavior of the form factors in the 
DSE approach is related to the behavior of the mass func¬ 
tion M(jr). At lower momenta, where the mass function is 
far from its perturbative or current quark value, the Pauli 
form factor is also large compared to its perturbative value 
and is falling more slowly than perturbation theory pre¬ 
dicts. (For reference, perturbative QCD predicts a Q -4 
power law falloff for Fi(Q 2 ) at large Q 2 , and a Q~ 6 falloff 
for F 2 (Q 2 ). ) Hence one can get a zero in Ge{Q 2 ), 

G e {Q 2 ) = F 1 {Q 2 )--^F 2 {Q 2 ), (51) 

and hence a falloff in the ratio Ge{Q 2 )/Gm(Q 2 )- 

For the newer DSE calculations reported by Segovia 
et al. [284], the zero in Ge is at Q 2 = 9.5 GeV 2 . If 
the mass function fell to its low perturbative value more 
quickly than it does, the quarks would behave more like 
free quarks, and the value of the Pauli or anomalous mag¬ 
netic moment, form factor would be small as well as quickly 
falling. In such a case, the zero of Ge{Q 2 ) would be pushed 
to higher values of Q 2 or possibly not occur at all 1284) . 

We may mention that models based on the Dyson- 
Sclrwinger equations, as well as many of the other models 
discussed, do extend to form factors for other hadronic 
reactions, such as the electromagnetic N —> A transi¬ 
tion [28411285] . The result for the ratio of the electric and 
magnetic transition form factors for this process, Rem-, 
turns out to be small in the DSE approach, even at mo¬ 
mentum transfers above 5 GeV 2 , in accord with experi¬ 
mental data. The perturbative QCD result, that Rem —> 
1, may well ensue, but only at momentum transfers well 
above those now experimentally accessible. 


4.3 Links between Deep-Inelastic Scattering and 
Nucleon Form Factors 


4.3.1 Perturbative QCD Inspired Models 

Perturbative QCD (pQCD) predicts the scaling behavior 
of the nucleon electromagnetic form factors are high Q 2 . 
The predictions were given by Brodsky and Farrar 12861 
and by Matveev, Muradyan, and Tavkheledze 


and pQCD not only predicts the scaling behavior of each 
amplitude, but also predicts that the leading amplitude 
is the one where the hadron helicity is maintained in the 
interaction. A photon of high Q 2 sees the nucleon as a 
set of three parallel moving quarks, one of which absorbs 
the high momentum photon, whose momentum is brought 
in from a sideways direction. In order to reconstitute the 
proton, the momentum must be shared out by two hard 
gluon exchanges, so that one also has three parallel moving 
quarks in the final state, as illustrated in the Fig. | 

The overall hard amplitude can be factored 
imil292ll293l^M| as a product of a hard scattering am¬ 
plitude that takes three parallel moving quarks into three 
parallel moving quarks, and two distribution amplitudes 
(DA) that specify how the longitudinal momentum of the 
nucleon is divided among the quarks. Each gluon carries 


N 



Fig. 29. Perturbative QCD picture for the nucleon electro¬ 
magnetic form factors. The highly virtual photon resolves the 
leading three-quark Fock states of the nucleon, described by 
a distribution amplitude. The large momentum is transferred 
between the quarks through two successive gluon exchanges 
(only one of several possible lowest-order diagrams is shown). 



a virtuality proportional to Q 2 (and there are also factors 
1 /Q from each of the internal quark propagators, factors 
of Q from each of the thoroughgoing quark lines, and a 
1 /Q involved in the definition of Fj), leading to a pQCD 
prediction that the helicity conserving Dirac form factor 
F\ will fall like 1/Q 4 , with possible logQ 2 factors, at high 
Q 2 . The Pauli form factor F 2 requires a helicity flip be¬ 
tween the final and initial nucleon, which in turn requires, 
thinking of the quarks as collinear, a helicity flip at the 
quark level, which is suppressed at high Q 2 . The result is 
a prediction that F 2 will fall like 1 /Q 6 at high Q 2 . Hence, 
Ge and Gm will both fall like l/Q 4 asymptotically. 
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Fig. 30. Test of the scaling behavior of the proton form factors. 
Upper panel: proton Dirac form factor multiplied by Q 4 . Lower 
panel: proton Pauli form factor multiplied by Q 6 . In both pan¬ 
els Ge p was obtained from the JLab fi p Ge p /Gm p data using 
the Kelly parametrization for Gm p 1551 . 


We can see how well pQCD predicts current electro¬ 
magnetic proton form factor data by examining Fig. 1301 
The figure shows data up to 10 GeV 2 for Q 4 Fi (upper 
panel) and Q 6 F 2 (lower panel). For Fi, it appears that 
the curve is flattening out, as pQCD would predict, and 
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indeed there is data for JFi up to 31 GeV 2 to corroborate 
this. However, for Q 6 F 2 , where there is no further data 
currently, the existing data does not match the simple 
pQCD expectation. 

The data show that F-ipj F\ p falls slower than 1 /Q 2 
with increasing Q 2 . Belitsky, Ji, and Yuan Esa investi¬ 
gated the assumption of quarks moving collinearly with 
the proton, which underlies the pQCD prediction. They 
have shown E25] that by including components in the nu¬ 
cleon liglit-cone wave functions with nonzero quark or¬ 
bital angular momentum projection, they obtain a behav¬ 
ior F 2 /F 1 —f In 2 (Q 2 / A 2 )/Q 2 at large Q 2 , with A a non- 
perturbative mass scale. (Refs. 129611297] also discuss using 
quark orbital angular momentum to get a ratio F 2 P /F\ p 
which drops slower than 1 /Q 2 with increasing Q 2 .) With 
A around 0.3 GeV [ 255] * the data for F 2 P /F\ p agree with 
such double-logarithmic enhancement, as seen in Fie. 1311 
where it may be noted that the higher Q 2 data was ob¬ 
tained after [295] was published. The arguments of [295] 
do still rely on pQCD, and it remains to be seen if still 
higher Q 2 data will continue to support this amended pre¬ 
diction. 

Hard scattering is calculated as if all three quarks 
are moving fast. An alternative is that one quark carries 
nearly all the nucleon’s momentum, and the other two 
quarks are soft. It is not necessary to transfer momen¬ 
tum to the soft quarks before reconstituting the proton. 
Nesterenko and Radyushkin [298] point out that the hard 
scattering mechanism requires the exchange of two gluons, 
each of which brings in a suppression factor a s / 7 r ~ 0 . 1 . 
One therefore see that the hard scattering mechanism for 
Ff could be numerically suppressed relative to the soft 
term, also called the Feynman mechanism; see also 12991 
19991 . 

Early on, Duncan and Mueller Ena. showed that the 
soft or Feynman process also gave a 1/Q 4 falloff, with log¬ 
arithmic corrections, for the Dirac form factor Fj. This 
has been taken up more recently by Kivel and Vander- 
haeghen [30211303] . who were able to show that also for 
the Feynman process, a type of factorization was possible, 
where the second step is given by a process independent 
kernel that transfers momentum among the initially all fi¬ 
nite momentum fraction quarks to make two of them soft, 
as in Fig. [32] They also considered the Pauli form factor 
F 2 , not with the same success in proposing a factorization 
theorem, but showing F 2 /Fi ~ 1 /Q 2 at high enough Q 2 
also for the soft process. 

We shall also mention work where the soft contribu¬ 
tion was evaluated within the light-cone sum rule (LCSR) 
approach of Braun et al. m- Using asymptotic distribu¬ 
tion amplitudes for the nucleon, the LCSR approach yields 
values of Gm p and Gmu which are within 20 % compatible 
with the data in the range Q 2 ~ 1 10 GeV 2 . The electric 
form factors however were found to be much more difficult 
to describe, with Geu overestimated, and Ge p /Gm p near 
constant when using an asymptotic nucleon distribution 
amplitude. Only when including twist-3 and twist-4 nu¬ 
cleon distribution amplitudes within a simple model, is a 
qualitative description of the electric proton and neutron 




Fig. 31. Upper panel: Data from the JLab plotted as 
Q 2 F 2 P /Fi p /In 2 (Q 2 /A 2 ) as proposed by Belitsky et al. 12951 . 
The solid curve is the fit to the data using Eqn. 44. Lower 
panel: the same but for the neutron, data from [32l33lHl2l 
mm - The solid curve is calculated from Kelly fit [55] . 


form factors obtained. Such higher twist components hint 
at the importance of quark angular momentum compo¬ 
nents in the nucleon wave function. 


4.3.2 Generalized Parton Distributions 

Generalized parton distributions (GPDs) represent an am¬ 
plitude for removing a quark from a nucleon and replac¬ 
ing it with another quark with a different momentum, 
and possibly also with different spin projection and fla¬ 
vor. These amplitudes can be measured in virtual Comp¬ 
ton scattering, r y*(qh) + N(p) —> 7 ( 9 ') + N(p'), or in meson 
electroproduction, e.g., 7 * + N —> p + N. The momenta 
are indicated, with q ' 2 = 0 , the virtuality Q 2 = — q 2 > 0 , 
and q will be the momentum transfer to the nucleon. 
A diagram for the virtual Compton process is shown in 
Fig. 1331 in a diagram where both photons interact with 
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the same quark. The upper part of the diagram is to be 
calculated perturbatively, and the lower part of the dia¬ 
gram is given by the GPD. For “deep” virtual Compton 
scattering, where Q\ Q 2 ,m%, it has been shown that 
the dominant contributions are given by diagrams like the 
one shown, and that one can separate or factor the per¬ 
turbative process specific stage of the interaction from the 
non-perturbative process independent part (see 3041 13051 
1505113071130811505] for reviews and references). 

Further notation is that P = (p + p')/ 2 is the average 
nucleon momentum, k is the average momentum of the 
quarks entering and leaving the nucleon, x is the light- 
front momentum fraction defined from x = k + /P + , and 
the asymmetry between the quark momenta is given by 
the skewedness £ = —q + / ( 2 P + ). In the high Q\ limit, one 
can show that £ is related to the Bjorken variable ig by 
2£ = x B /(1 - x B / 2), where x B = Q 2 h /{2p ■ q h ). 

GPDs were introduced by Ji [3101 andjpy Radyushikin 
m- Formally, in the notation of Ji 131 (ifl and in a frame 
where P and qh are collinear with P in the positive z- 


1 A gauge link Pexp(igr f dx M A M ), ensuring color gauge in¬ 
variance, is tacit. 



Fig. 32. Soft contributions to the nucleon form factor. The 
kernels J connect three quarks with finite momentum to a 
configuration with two soft quarks and a hard quark interacting 
with the photon. In the simplest case, each kernel J would 
contain two gluon exchange. 



Fig. 33. The “handbag” diagram for the nucleon DVCS pro¬ 
cess. Provided the virtuality of the initial photon (with mo¬ 
mentum qh) is sufficiently large, the QCD factorization theo¬ 
rem allows to express the total amplitude as the convolution of 
a Compton process at the quark level and a non-perturbative 
amplitude parameterized in terms of generalized parton dis¬ 
tributions (lower blob). The diagram with the photon lines 
crossed is also understood. 


direction, one obtains the GPDs from, 



x (N(p'M q (-y/2) 7+ *l>q{y/2)\N(p))\ y+ ^ y±=0 
= H q (x,tQ 2 ) N{p')j + N{p) 

+ E q (x, £, Q 2 ) N(p) ^^ v q v Nip), (52) 

where if> q is the quark field for flavor q and N is the nucleon 
Dirac spinor. The matrix element is non-perturbative and 
is given in terms of two functions H q and E q for each 
flavor q. There is a similar matrix element with operator 
V’g7 + 75V , g and two further, polarized, GPDs H q and E q . 

The notation is such that positive momentum fractions 
correspond to quarks and negative ones to antiquarks. 
Hence in a; > £, the fermions leaving and entering the 
nucleon are both quarks. Also possible is that both are 
antiquarks (x < — £), or that there is a quark-antiquark 
pair. 

GPDs enter this review because of their relation to 
form factors. An integral over x on the LHS of Eq. (1551) 
forces both quark fields to be at the same point, as in the 
matrix elements of the electromagnetic current. One can 
show that the matrix elements have support for — 1 < x < 
1 , and that 

j + \xH q (x : t 1 Q 2 )=F?(Q 2 ), 

I* dxE q (x,Z,Q 2 ) = F q (Q 2 ), (53) 

where the nucleon form factor are given in terms of the 
quark flavor form factors F® in the expected ways, 

p. — - p u _ i p d _ i p s 
ip — 3 ± 1 3 ± 1 3 1 ’ 

F in = -±F? + lF l d -±F°, (54) 

where * = 1,2 and F™’ d are specifically for the proton. 

These relations allow us, if we have complete mea¬ 
surements or good models (the latter is more the case 
at present) for the GPDs, to obtain the electromagnetic 
form factors from them. Alternatively, the measured form 
factors can be used as constraints upon GPD models. 

An example of model GPDs is the modified Regge pa¬ 
rameterization for H and E that was proposed by Guidal 
et al. 3121 . 

H q (x,0,Q 2 ) = q v {x)x a ' C 1 -*)® 2 , 

E q (x,0,Q 2 ) = (1 - x)^ q v {x)x a ' (1 ~ x)q2 , (55) 

depending on 3 parameters. The Regge slope a' is deter¬ 
mined from the Dirac radius, and two parameters q u and 
rj d , entering the GPD E, ensure that the x ~ 1 limit of E q 
has extra powers of 1 — x compared to that of H q . This 
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Q 2 (GeV 2 ) 


Fig. 34. GPD calculation of Gm p relative to the dipole Gd 
and of Gep/Gmp in the upper panels, and the correspond¬ 
ing plots for the neutron in the lower panels, using GPD’s 
from Ref. [3121 . Data for Gm p are from [12] (open squares), 
ITT , (open circles), |T5] (blue solid stars), [TB] (green open 
stars), [20 ] (red solid circles), [54] (red solid squares), accord¬ 
ing to the recent re-analysis of Ref. [56]. Data for the ratio 
Gep/Gmp are from [02] (blue open triangles), [5] (red solid 
squares), :8] (blue solid circles), and llOOl (green solid tri¬ 
angles). The data for Gmu are from H3lj (red solid circles), 
1132] (red solid squares), [76] (open triangles), [77] (green open 
stars), [72] (open squares), 1313] (solid triangles), and [314] 
(blue solid stars). The data for Geti are from double polariza¬ 
tion experiments at MAMI |4lll42llTTTllll8llll9l (red solid cir¬ 
cles), NIKHEF [TO] (green solid triangle), and JLab |31II112I 
132'] (blue solid squares). 


results in a proton helicity flip form factors F<i which has 
a faster power fall-off at large Q 2 than F\, as observed 
experimentally. 

In Fig. [34] the proton and neutron Sachs electric and 
magnetic form factors are shown. One observes that the 
3-paranreter modified Regge model gives a rather good 
overall description of the available form factors data for 
both proton and neutron in the whole Q 2 range, using as 
value for the Regge trajectory a' = 1.105 GeV -2 , and the 
following values for the coefficients governing the x —> 1 
behavior of the .E-type GPDs: rj u = 1.713 and rj d = 0.566. 
Note that a value rj q = 2 corresponds to a 1/Q 2 asymp¬ 
totic behavior of the ratio F%/F^ at large Q 2 . The modi¬ 
fied Regge GPD parameterization allows one to accurately 
describe the decreasing ratio of Ge p /Gmp with increasing 
Q 2 , and also leads to a zero for Ge p at a momentum trans¬ 
fer of Q 2 ~ 8 GeV 2 . 


4.4 Lattice QCD Calculations of Nucleon Form Factors 

Strictly speaking, lattice calculations of nucleon form fac¬ 
tors are currently available only for the isovector form fac¬ 
tors. 

Isoscalar form factors require calculations of discon¬ 
nected diagrams, which are diagrams with quark loops not 
connected to the quark lines emanating from or ending on 
the lattice nucleon source or sink. There are gluons that 
attach the quark loops to the valence quarks, but these are 
not indicated in lattice diagrams, hence the phrase “dis¬ 
connected.” Contributions from the disconnected loops re¬ 
quire computer time intensive calculations, and remain 
undone. However, the disconnected diagrams contribute 
equally to proton and neutron, so the isovector case can 
be considered without them. 

A review including lattice form factor results up to 
2010 is available in [3T5] . and newer lattice form factor 
results are reported in 131611317(1318] . 

The new calculations reported in Green et al. [318] 
have pion masses from 373 MeV down to close to phys¬ 
ical 149 MeV. The latter also strove to reduce contami¬ 
nation from excited nucleons. They analyze their lattice 
data using three methods which they call the standard 
ratio method, the summation method, and the general¬ 
ized pencil-of-function method (GPoF), with varying out¬ 
comes. The best results, judged by comparison to data as 
represented by one of the standard fits I3l9j . come from 
the summation method. Here agreement with experimen¬ 
tal data is good for both Ge v and Gmv in the region 
considered, which is Q 2 from scattering threshold up to 
about 0.5 GeV 2 , with uncertainty limits about 20% at Q 2 
of 0.4 GeV 2 . 

The works of Alexandrou et al. [ 3161 and Bhattacharya 
et al. m have pion masses in the 213-373 MeV range, 
and quote results for somewhat higher Q 2 . For Q 2 above 
about 0.6 GeV 2 , their isovector form factors results tend 
to be 50% or so above the data for Ge v (or Fi v ), with 
uncertainties indicated at about 10%. For Gmv (or F 2 V ) 
the results are closer to data . The authors of these works 
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do point out that the lattice treatments with these pion 
masses are all consistent with each other. 

One may specifically focus on nucleon radii calculated 
from lattice gauge theory. In the future, it may be possi¬ 
ble and desirable to calculate using a dedicated correlator 
which gives directly the slope of the form factor at zero 
momentum transfer. Finding such correlators by taking 
derivatives of known correlators is suggested and stud¬ 
ied 1320J for lattice calculations of form factors at points 
where the Lorentz factors they multiply go to zero. Ap¬ 
plications in nreni are to form factors for semi-leptonic 
scalar meson decay, and to hadronic vacuum polarization 
corrections to the muon (g — 2 ). 

At present, lattice calculations of nucleon radii pro¬ 
ceed by calculating the form factor at several non-zero 
Q 2 , fitting to a suitable form, typically a dipole form, and 
finding the radius by extrapolating to zero Q 2 . Truly com¬ 
plete results are available only for the isovector nucleon. 
Ref. [3181 presents a plot of radius results for lattice calcu¬ 
lations at various pion masses. They use the Dirac radius, 
obtained from the slope of F\ v , rather than the charge ra¬ 
dius obtained from Ge v , but these are related by, using 
the proton as an example, 

= ( 56 ) 

where k v is the isovector anomalous nucleon magnetic mo¬ 
ment. Hence, given the great accuracy of the magnetic 
moment measurements, one knows the Dirac radii to the 
same accuracy as the charge radii. 

The great interest is to obtain sufficient accuracy from 
the lattice results to be able to adjudicate between the 
electron and muon measured values of the isovector charge 
or Dirac radii. The electron measured isovector radius is 
straightforward to look up, the muon measured value of 
the Dirac or charge radius is for now a defined quantity ob¬ 
tained by using the electron value for the neutron radius- 
squared. Using the summation method, Ref. [3T8j obtains, 
by extrapolation to the physical pion mass, a value of 
the isovector Dirac radius between the muonic and elec¬ 
tronic results, with uncertainties that accommodate both 
at about the one standard deviation level. However, using 
the GPoF or ratio method gives a smaller (rf) v , on the 
order of 2/3 the value from the summation method. 

One may say there is opportunity for further work. An 
uncertainty of 1 % or less for the proton alone is needed for 
a lattice calculation to impact the proton radius puzzle. 
Two extrapolations are needed to obtain the charge ra¬ 
dius for a physical nucleon. On is in the lattice pion mass, 
commented upon above. The other is in Q 2 . Currently on 
the experimental side, the lowest Q 2 from scattering data 
is about 0.004 GeV 2 , and there are experiments planned 
or running to reduce this number. Further there are dis¬ 
cussions, alluded to in Sec. 14.1.11 regarding the best fit 
forms to use for the extrapolation. Currently on the lat¬ 
tice, represented by Ref. EH], the lowest Q 2 is about 0.04 
GeV 2 , and the fit to the lattice data is only done using a 
simple dipole form. 



Fig. 35. The isovector electric form factor G v Ep = Ge p -Gb„ 
obtained from fits to the experimental form factor data shown 
in Figs. IKil and ED 


And of course one wants the isoscalar as well as the 
isovector form factors. However, a challenge involving just 
the isovector form factors, albeit at a higher Q 2 than lat¬ 
tice form factor results currently display, is to obtain the 
zero in the isovector Ge v that is visible in the fits to the 
experimental data shown in Fig. 1351 


5 Outlook 

The experimental and theoretical status of the nucleon 
form factors were reviewed extensively in the 15 years 
following publication of the results of the first recoil po¬ 
larization experiment at Jefferson Lab [7] for the proton. 
Chronologically these reviews include Gao m, Hyde- 
Wright and de Jager 13251 . Perdrisat et al. [3231, Arrington 
et al. [324], Cloet et al. |123] . Arrington et al. |325] . Per¬ 
drisat and Punjab i EH and S. Pacetti et al. [327| . The 
completion of the GEn(l) and GEp(3) experiments, which 
reached a maximum Q 2 of 3.4 and 8.4 GeV 2 respectively, 
has brought the field into previously unexplored regions 
of four-momentum transfer squared, and correspondingly, 
generated a burst of theoretical investigations along old 
and new paths. 

The proton form factors were originally introduced in 
the approximation of non-relativistic scattering, as the 
three-dimensional Fourier transform of the charge den¬ 
sity [251l328| . However the proton recoil implies that the 
electron is interacting with a moving charge distribution. 
Already for Q 2 =0.25 GeV 2 , the recoil proton relativistic 
boost factor 7 is 1.1, corresponding to v/c = 0.42. The 
argument that form factors are Fourier transforms of nu¬ 
cleon density in the Breit frame had to be abandoned, 
as this frame’s velocity in the Lab frame is significantly 
different for every Q 2 . 

Changes in our view of the structure of the proton are 
many. For example, the proton in its ground state is not 
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Fig. 36. Upper panel: Five theoretical model predictions of the 
proton form factor ratio [i p Ge p IGm p shown together with the 
data of JLab recoil polarization experiments, symbols are same 
as in Figure 10. Lower panel: The same five theoretical model 
predictions for the neutron form factor ratio fi n GE P /GMn, 
data are same as in Figure 11. The theoretical curves are from 
the VMD model of Lomon |208[ , solid (black) curve; from the 
GPD’s of Guidal et al 1312| . medium dashes (magenta), the 
covariant spectator model of Gross et al. 1239j as short dash- 
dot (green) Dyson-Schwinger equation solutions of Cloet et al. 
12821 . long dash-dots (cyan), and the quark-diquark model of 
Cloet and and Miller [260] . medium dash-dots (blue). The sym¬ 
bols for the data are explained in Fig. 1111 


necessarily spherically symmetric, but can show a typical 
multipole shape, when referred to the spin direction of one 
of its quarks (constituents) relative to the nucleon spin ori¬ 
entation 13231 . Also, the wave front or infinite momentum 
frame charge and magnetization densities are invariant, 
two-dimensional transverse distributions which are dras¬ 
tically different from the non-relativistic ones [32911330] . 

A selection of model predictions for the form factor 
ratios ii n GEn/GMn and h p Ge p IGm p is shown together 
with the data obtained in double polarization experiments 
at JLab and in other laboratories in Fig. [351 These two 
figures emphasize the importance of future experiments 
which will establish whether either ratio does, or does not 
cross zero near Q 2 =10 GeV 2 . Whereas for h p Ge p /Gm p , 
all model predictions discussed in this review anticipate a 
zero crossing somewhere above 9 GeV 2 , for n n GEn/GMn 
only the calculation based on the Dyson-Schwinger equa¬ 
tions predict such a zero crossing in the 10 GeV 2 region 
of Q 2 . 

A recent development has been the calculation of the 
flavor separated form factors of the ’’dressed” quarks from 
simple linear relations between the nucleon form factors, 
assuming charge symmetry applied to the data available. 
The dressed up and down quarks have significantly differ¬ 
ent form factors |123lll91[l331H332ll2821ll92j . Nucleon form 
factors determine the parameters of the valence quark 
GPDs; these can be used to obtain corresponding valence 
quark densities 13331 . They can be compared with the 
GPDs obtained from real and virtual Compton scatter¬ 
ing. 

The doubling of the energy of the Jefferson Lab accel¬ 
erator to 12 GeV will lead to a much enhanced program 
of experiments investigating the structure of the nucleon. 
An experiment will use the existing High Resolution Spec¬ 
trometers in Hall A at Jefferson Lab to measure Gm p with 
greatly improved error bars up to 14 GeV 2 [ 3331 . An new, 
versatile Super BigBite Spectrometer (SBS), consisting of 
a simple dipole magnet and associated detectors, is be¬ 
ing built for three form factor experiments in Hall A. One 
SBS experiment will measure Ge p /Gm p up to Q 2 = 12 
GeV 2 using the recoil polarization technique 1334] . Sitting 
behind the SBS dipole magnet will be a recoil polarimeter 
which will have two analyzers with multiple GEM cham¬ 
bers used for incoming and scattered track determination. 
Another SBS experiment will extract Geu/Gmu up to 
Q 2 = 10 GeV 2 from beam-target asymmetry measure¬ 
ments using an upgraded polarized 3 He target [3351 . The 
experiment will detect the scattered electrons in the Big- 
Bite spectrometer and the scattered neutron in a large 
solid angle hadron calorimeter sitting behind the SBS mag¬ 
net. In Hall C, measurements of the neutron recoil polar¬ 
ization in quasi-free electron deuteron scattering will be 
done to extract Geu/Gmu up to Q 2 = 6.9 GeV 2 336 . 
With the new CLAS12 spectrometer in Hall B at Jef¬ 
ferson Lab , the measurement of Gmu will be done to 
Q 2 = 14 GeV 2 [ 3371 . A third SBS experiment will also 
measure Gmu to Q 2 = 14 GeV 2 [3381 . The proposed error 
bars for all these experiments are shown in Fig. 1371 
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Fig. 37. The projected error bars for the approved nucleon form factors experiments at Jefferson Lab in the 12 GeV era. For the 
Hall A SBS experiment E12-07-108 13341 . the anticipated error bars on the ratio h p Ge p /Gm p are shown as filled circles (cyan). 
The anticipated error bars on the ratio Geti/Gm™ are shown as filled squares (magenta) for the Hall A SBS experiment E12-09- 
016 [335J and as filled diamonds (cyan) for the Hall C experiment E12-11-009 [336] , The anticipated error bars for Gm p /h p Gd 
from the Hall A experiment E12-09-019 13331 are shown with square symbols (magenta). Finally the ratio Gmti/HuGd will be 
measured in two experiments: E12-09-019 in Hall A [3381 and E12-07-104 in Hall B 13371 . The expected error bars are shown as 
empty circles (magenta) and filled stars (cyan), respectively. 


One of the most stringent constraints that nucleon 
elastic form factor data at large Q 2 can provide, relates 
to the issue of the various contributions from quarks, glu¬ 
ons, and orbital angular momentum to the total angular 
momentum of the nucleon. The elastic form factors also 
provide a powerful check of lattice QCD. The lattice calcu¬ 
lations of form factors are making impressive progress, and 
the comparison of these results with experimental mea¬ 
surements will be extremely important. There is an indi¬ 
cation from the results of GEp(3) experiment that we may 
be entering the range of momentum transfers where the 
pQCD prediction is vindicated. Yet a continuation of the 
fast decrease of the ratio toward negative values cannot 
be excluded. Great progress in the theoretical description 
of the structure of the nucleons can be expected. 
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